Projects  and  Publications 

of  the 

NATIONAL  APPLIED  MATHEMATICS  LABORATORIES 


A QUARTERLY  REPORT 
July  through  September  1949 


NATIONAL  APPLIED  MATHEMATICS  LABORATORIES 

of  the 

NATIONAL  BUREAU  OF  STANDARDS 


national  applied  mathematics  laboratories 

July  1 through  September  30,  1U49 


ADMINISTRATIVE  OFFICE 

John  II.  Curtiss,  Ph.D.,  Chief 
Edward.  W.  Cannon,  Ph.D.,  Assistant  Chief 
Mrytle  R.  Kellington,  M.A.,  Technical  Aid 
Luis  0-  Rodriguez,  M.A.,  Chief  Clerk 
John  3.  Tallerico,  B.C.S.,  Assistant  Chief  Clerk 
Jacqueline  Y.  Kartell,  Secretary 
Dora  P.  Cornwell,  Secretary 
Esther  McCraw,  Secretary 
Pauline  F.  Peterson,  Secretary 


INSTITUTE  FOR  NUMERICAL  ANALYSIS 
Los  Angeles,  California 


COMPUTATION  LABORATORY 


John  H.  Curtiss,  Ph.D Acting  Chief  ( July-August ) 

J.  Earkley  Rosser,  Ph.  D.  ...  Director  of  Research  (September-  ) 
Albert  S.  Cahn,  Jr.,  M. S Asst,  to  Director  of  Research 

Research  Staff 


Forman  S.  Acton,  Ph.  D.  .....  . 

Lars  V.  Ahlfors,  Ph.  D 

Edwin  F.  Beckenbach,  Ph.D. . . . 

Monroe  D.  Donsker,  Ph.D 

William  Feller,  Ph.D 

George  E.  Forsythe,  Ph.D 

Samuel  Herrick,  Ph.D 

Magnus  R.  Hestenes,  Ph.D 

Mark  Kac,  Ph.D 

William  Karush,  Ph.D 

Cornelius  Lanczos,  Ph.D 

Alexander  M.  Ostrowski,  Ph.D. 

Otto  Szasz,  Ph.D 

Stephen  E.  Warschawski , Ph.D. 

Wolfgang  R.  Wasow,  Ph.D 

Graduate  Fellows: 

Robert  J.  Diamond,  B.S. 

Robert  C.  Douthitt,  M.A. 
Ernest  S.  Elyash,  B.S. 

George  Gourrich,  B.S. 

Harold  Gruen,  M.A. 

Lloyd  K.  Jackson,  M.A. 


Mathematician 

...  Consultant  (Mathematician) 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathe.nat  ician 

. . . Consultant  (Mathematician) 

Mathematician 

Mathematician 

Mathematician 

Raymond  P.  Peterson,  Jr.,  M.A. 

Marvin  L.  Stein,  M.A. 

James  G.  C.  Templeton,  M.A. 

Hans  F.  Weinberger,  M. S. 

Robert  Sehnert,  B.S. 


Machine  Development  Unit 


Harry  D.  Huskey,  Ph.D 

Biagio  F.  Ambrosio,  B.S 

Edward  L.  Lacey,  M. S 

Roselyn  Siegel  Lipkis,  B.A 

David  F.  Rutland,  M. S 

Laboratory  Electronic  Mechanics : 
Robert  V.  Barnett 
Charles  M.  Brauer 
Arnold  Dolmatz 
Blanche  C.  Eidem 
Sidney  S.  Green,  B.S. 

Computation  Unit 


Mathematician 

. Electronic  Engineer 
. Electronic  Engineer 

Mathematician 

Electronic  Scientist 

Harry  T.  Larson,  B.S. 

John  D.  Mack 
John  L.  Newberger 
Geraldine  Orr 
James  W.  Walsh 


Gertrude  Blanch,  Ph.D.... 
Alfredo  Banos,  Jr.,  Ph.D. 
David  S.  Saxon,  Ph.D 
Louis  L.  Bailin,  Ph.D 
Marvin  Howard 
Everett  C.  Yowell, 
Computers: 

Helen  C.  Arens,  M.A. 
Patricia  Burton  Bremer, 
Leola  Cutler,  B.A. 

Lillian  Forthal,  B.A. 
Fannie  M.  Gordon,  M.B.A. 
Elizabeth  F.  Harding,  B.A. 
Donald  E.  Henley 
Ruth  B.  Morgan,  B.A. 

Gerald  W.  Kimble,  A.B. 

Eileen  Feyder 

M.  Winifred  Smith 


Mathemat ician 

Consultant  (Mathematician) 
Consultant  (Mathematician) 
Mathematician 
Mathematician 
Mathematician 

Nancy  Robbins  Mann,  M.A. 
Shirley  L.  Marks,  M.A. 
William  0.  Paine,  Jr. 
John  A.  Postley,  B.A. 
Everett  A.  Rea,  B.A. 
Albert  H.  Rosenthal 
Herbert  Snow,  B.A. 
Rose  Tishman.  B.A. 
Mary  J.  Tudor,  B.A. 

Key  Punch  Operator 

Key  Punch  Operator 


M.S 

Ph.D 

B.A. 


Wilbur  W.  Bolton,  Jr. 
M.  Evelyn  Michaud.... 
F.  Patricia  K.  Peed.. 
Ellma  Frantz,  B.S.  . . 
Velma  R.  Huskey,  B.A. 

Secretaries ; 

Elsie  E.  Aho 
Irene  M.  Everett 
Mildred  S.  Goldberg 
Margaret  W.  Gould 
Eleanore  J.  Harris 


Administrative  Unit 

B.A Adm.  Officer  ( July-September ) 

..  Administrative  Officer  (September-  ) 

Librarian 

Library  Assistant 

Physical  Science  Editor 

Dorothy  M.  Hibbard 
Elsie  L.  Kusman 
Dora  Kaplan 
Sata  Ohashi 
Gertrude  Z.  Reider 


Franz  L.  Alt,  Ph.D 

Arnold  N.  Lowan,  Ph.D 

John  Todd 

Jack  Belzer,  B.A 

Natalie  Coplan,  B.A 

Wanda  P.  Gordon,  B.A 

Benjamin  F.  Handy,  Jr.,  M.S. 

Joseph  B.  Jordan,  B.A 

Henry  B.  Juenemann,  B.A. . . . 

Joseph  H.  Levin,  Ph.D 

Michael  S.  Montalbano,  B.A. . 
Malcolm  W.  Oliphant,  M.A. . . 

Irene  A.  Stegun,  M.A 

Milton  Stein  

Ruth  Zucker,  B.A 

Computers: 

John  A.  Astin 
Ruth  E.  Capuano 
Bernard  C.  Dove 
Elizabeth  F.  Dodefroy 
Viola  D.  Hovsepian 
Lionel  Levinson 
David  S.  Liepman 

Oneida  L.  Baylor 

Florence  C.  Pettepit 

Lillian  Sloane 


...  Asst,  and  Acting  Chief 
Consultant  (Mathematician) 
Consultant  (Mathematician) 

Mathemat ician 

Mathematician 

Mathematician 

...  Gen'l  Phys ’ 1 Scientist 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mary  Orr 
B.  Stanley  Prusch 
Andress  0.  Ridgway 
Ruth  C.  Shepherd 
Jack  Singer 
Bertha  H.  Walter 
Marjorie  0.  White 

Card  Punch  Operator 

Secretary 

Secretary 


STATISTICAL  ENGINEERING  LABORATORY 


Churchill  Eisenhart,  Ph.D. 

W.  J.  Youden,  Ph.D 

Albert  H.  Bowker,  Ph.D..., 
Joseph  M.  Cameron,  M.S. ... 

Lola  S.  Deming,  M.A 

Gerald  J.  Liebcrman,  M.A.. 

Julius  Lieblein,  M.A 

Celia  S.  Martin,  B.A 

Vivian  M.  Frye,  B.A 

Helen  M.  Herbert 


Chief 

Assistant  Chief 

Consultant  (Mathematician) 

Mathematician 

Mathematician 

Math’l  Statistician 

Mathemat ician 

Mathematician 

Secretary 

Secretary 


MACHINE  DEVELOPMENT  LABORATORY 


Edward  W,  Cannon,  Ph.D. 
Samuel  Lubkin,  Ph.D.... 

Ira  C.  Diehm 

Florence  K.  Koons,  B.A. 
Ethel  C.  Marden,  B.A. . . 
Edith  T.  Norris,  B.A... 

Ida  Rhodes,  M.A 

Otto  T.  Steiner,  B.A... 
Susan  A.  Moon 


Acting  Chief 

Electronic  Scientist 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Mathematician 

Secretary 


U.  S.  DEPARTMENT  OF  COMMERCE 
Charles  Sawyer,  Secretary 


NATIONAL  BUREAU  OF  STANDARDS 
E.  U.  Condon,  Director 


Projects  and  Publications 

of  the 

NATIONAL  APPLIED  MATHEMATICS  LABORATORIES 

July  through  September  1949 


Contents 


Page 


Index iii 

Status  of  Projects  as  of  September  30,  1949 1 

Institute  for  Numerical  Analysis  (NBS  Section  11.1).  . 1 

Computation  Laboratory  (NBS  Section  11.2)  29 


Statistical  Engineering  Laboratory  (NBS  Section  11.3)  . 43 

Machine  Development  Laboratory  (NBS  Section  11.4)  ...  54 


Lectures  and  Symposia.  . . 59 

Publication  Activities  6l 


National  Applied  Mathematics  Laboratories  of  the  National  Bureau  of  Standards 


Preface 


This  is  a report  on  the  activities  of 
Division  11  of  the  National  Bureau  of  Standards 
for  the  period  from  July  1,  1949  through 
September  30,  1949. 

Division  11  is  known  as  the  National  Applied 
Mathematics  Laboratories.  It  is  the  mission  of  the 
Laboratories  to  perform  research  and  to  provide 
services  in  various  quantitative  branches  of  mathe- 
matics, placing  special  emphasis  on  the  development 
and  exploitation  of  high-speed  numerical  analysis 
and  modern  statistical  methodology.  The  Labora- 
tories maintain  an  expert  computing  service  of 
large  capacity,  and  provide  consulting  services  in 
classical  applied  mathematics  and  in  mathematical 
statistics.  These  services  are  available  primarily 
to  other  federal  agencies,  but  under  certain  cir- 
cumstances it  is  possible  to  perform  work  for 
industrial  laboratories  and  universities. 

Inquiries  concerning  the  availability  of  the 
services  of  the  National  Applied  Mathematics  Labora- 
tories, or  concerning  further  details  of  any  of  the 
projects  described  in  this  report,  should  be  ad- 
dressed to  the  National  Applied  Mathematics  Labora- 
tories, 415  South  Building,  National  Bureau  of 
Standards,  Washington  25,  D.  C. 


Director 

National  Bureau  of  Standards 
November  15,  1949 


Index  of  Active  Research  and  Development  Projects 


Note:  This  index  is  not  intended  to  cover  the  numerous  special  problem 

solutions,  statistical  analyses,  and  other  ad  hoc  services  to  Government 
agencies  which  form  an  important  part  of  the  work  of  the  National  Applied 
Mathematics  Laboratories.  These  services  are,  however,  fully  represented 
in  the  body  of  the  report. 

A.  Research:  Pure  Mathematics 

Miscellaneous  studies  in  pure  mathematics  16 

Orthogonal  polynomials.  Determinants  of  . 11 

B.  Research:  Numerical  Analysis 

Characteristic  functions  and  characteristic  values  of  arbitrary  linear 

differential  operators.  Approximate  determination  of  the 9 

Characteristic  values  of  matrices.  Determination  of  1 
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Critical  points,  Determination  of  ....  . 12 
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Calculation  of 13 
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Solution  of  sets  of  simultaneous  algebraic  equations  and  techniques 
for  the  inversion  and  iteration  of  matrices 2 

C.  Research:  Applied  Mathematics,  Physics,  and  Astronomy. 
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Orbits,  Theoretical  and  practical  problems  in  the  application  of 
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value  of  a simple  quantum  model 27 

Theoretical  amplitude  distortion  in  a condenser  microphone  circuit.  4 

D.  Mathematical  Statistics 

Interpolated  medians.  Statistical  procedures  for.  ....  44 

Mean  deviation,  standard  deviation,  and  range  as  estimators  of 
scale  parameters  (measures  of  dispersion)  of  probability  distri- 
butions, The 43 

Probability  theory  of  a certain  double  branching  system  51 

Samples  of  three  observations.  Statistical  properties  of 43 

Stochastic  processes.  Elementary  theory  of 44 

2x2  tables.  Statistical  tests  of  significance  for,  when  the 
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Status  of  Projects 

September  30,  1949 


I.  Institute  for  Numerical  Analysis 
(Section  11.1) 


1 • Research  In  Numerical  Analysis 


Key  to  letter  symbols  in  project  numbers: 


AE 

AM 

CM 

ODE 

PD 

PM 

SF 


Algebraic  equations 
Analytical  applied  mathematics 
Numerical  methods  in  conformal  mapping 
Ordinary  differential  equations 
Partial  differential  equations 
Probabilistic  methods  in  numerical  analysis 
Studies  of  special  functions 


DETERMINATION  OF  CHARACTERISTIC  VALUES  OF  MATRICES 
Project  11.1/1-49-AE1 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Manager:  C.  Lanczos 

Objective:  To  develop  a practical  and ’ economical  method  of  evaluat- 

ing the  characteristic  values  of  arbitrary  complex  matrices. 

Background:  The  calculation  of  characteristic  values  of  matrices 

is  of  importance  in  many  problems  in  physics,  engineering,  and  applied 
mathematics.  In  particular,  characteristic  values  play  an  important  role 
in  flutter  and  other  vibration  analysis  problems,  and  it  was  the  inade- 
quacy of  methods  currently  in  use  for  such  problems  which  furnished  the 
motivation  for  the  present  project. 

Comments:  The  principle  underlying  the  present  investigation  con- 

sists in  iterating  a trial  vector  to  set  up  a system  of  linear  equations 
having  a recurrent  matrix,  and  then  designing  an  algorithm  for  the  solution 
of  the  recurrent  system.  The  solution  of  the  circular  system  yields  the 
coefficients  of  the  characteristic  equation.  This  method  is  analogous 
to  the  method  on  which  project  11 . 1/1-49-ODE4  is  based. 

Status:  TERMINATED.  It  was  decided  that  the  scope  of  the  objective 

should  be  widened,  so  project  11.1/1-50-3  was  set  up  to  replace  this 
project.  See  current  status  report  for  project  11.1/1-50-3. 


Authorized  2/15/49 
Terminated  9/30/49 


1 
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Status  of  Projects 


SOLUTION  OF  SETS  OF  SIMULTANEOUS  ALGEBRAIC  EQUATIONS  AND 
TECHNIQUES  FOR  THE  INVERSION  AND  ITERATION  OF  MATRICES 
Project  11 . 1/1-49-AE2 

Origin:  NBS  * Authorized  2/15/49 

Sponsor:  Office  of  Naval  Research,  USN  Revised  9/15/49 

Managers:  J.  B.  Rosser,  C.  Lanczos,  Gertrude  Blanch, 

G.  E.  Forsythe 

Objective:  To  develop  practical  and  economical  methods  for  solving 

linear  and  non-linear  simultaneous  equations  in  many  unknowns,  to  be 
used  in  connection  with  various  types  of  computing  equipment,  and  to 
study  related  problems  of  matrix  inversion  and  iteration. 

Background:  Sets  of  simultaneous  algebraic  equations  and  related 

matrix  problems  arise  repeatedly  in  all  branches  of  applied  mathematics. 

Status:  CONTINUED.  Dr.  Lanczos'  results  are  Incorporated  in 

his  comprehensive  work  on  the  eigenvalue  problem;  see  project  11.1/1-50-5 
Dr.  Blanch  studied  various  techniques  for  "positioning"  the  elements  of 
a matrix.  She  will  pursue  this  matter  further  during  the  remainder  of 
the  year . 

In  general, 'work  on  this  project  was  subordinated  to  other 
interests  of  the  staff.  However,  at  the  end  of  the  quarter,  a conference 
on  the  project  was  held  by  Dr.  Rosser  and  greatly  increased  activity 
Is  expected  in  the  future. 


SOLUTION  OF  EQUATIONS  IN  A SINGLE  UNKNOWN 
Project  11 .1/1-49-AE3 

Origin:  NBS  Authorized  4/1/49 

Sponsor:  Office  of  Naval  Research,  USN  Revised  9/15/49 

Manager:  A.  M.  Ostrowski 

Objective:  To  find  new  and  better  ways  of  dealing  with  algebraic 

and  transcendental  equations  In  a single  unknown. 

Background:  Such  equations  arise  repeatedly  in  all  branches  of 

applied  mathematics. 

Status:  CONTINUED.  Dr.  Ostrowski  prepared  a note  giving  the 

results  obtained  in  a study  of  Vincent's  method  for  separating  the  real 
roots  of  an  algebraic  equation  (Publication  (2)1.  A summary  of  the 
results  is  included  in  the  previous  quarter's  sbatus  report.  The  old 
and  rare  textbook  mentioned  was  obtained,  but  did  not  contribute  any- 
thing new  to  the  project  manager's  earlier  knowledge  of  the  history  of 
the  method. 

Dr.  Ostrowski  continued  to  work  on  the  improvement  of  New- 
ton's method  for  calculating  the  roots  of  a polynomial. 

Publications:  (l)  "On  two  problems  in  abstract  algebra  connected 

with  Horner's  rule, "by  A.  M.  Ostrowski;  submitted  to  the  American 
Mathematical  Monthly.  (2)  "Note  on  Vincent's  theorem,"  by  A.  M. Ostrowski 
submitted  to  the  Annals  of  Mathematics. 


Status  of  Projects 


3 


SOLUTION  OF  THE  TELEGRAPHER'S  EQUATION  FOR  INITIAL 
CONDITIONS  GIVEN  ON  ONLY  ONE  CHARACTERISTIC 
Project  11.1/1-49-AM2 


Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Managers:  G.  E.  Forsythe,  Gertrude  Blanch 


Authorized  11/1/48 
Terminated  9/30/49 


Objective:  (l)  To  solve  the  telegrapher's  equation 


,2 

<)  v + 1 

ax  at  ¥ 


V 


o 


for  t > 0 with  the  sole  boundary  condition  that  v is  known  for  all  x when 
t = 0.  Here  x is  essentially  terrestrial  longitude  and  t is  time.  (2)  To 
investigate  the  solution  of  a related  difference  equation.  (3)  To  present 
the  solution  in  such  a form  that  meteorologists  can  see  its  consequences 
and  try  it  out  in  forecasting. 

Background:  This  boundary-value  problem  arises  in  Rossby's  linear- 

ized treatment  of  the  equation  of  non -divergent  flow.  The  harmonic  solu- 
tions have  been  discussed  many  times  and  it  is  considered  important  to  see 
what  happens  with  arbitrary  initial  conditions.  In  general,  the  problem 
has  no  unique  solution  for  - «s  <.  x <.  . The  fact  that  initial  conditions 

are  given  on  only  one  characteristic  makes  it  a most  unusual  problem  for 
a hyperbolic  partial  differential  equation.  In  the  present  treatment  the 
solution  is  rendered  unique  by  using  the  roundness  of  the  earth,  i.e., 
by  making  the  solution  periodic  in  x. 

Status:  TERMINATED.  Although  Dr.  Forsythe  has  certain  unpublished 

results  on  (2)  in  the  Objective,  there  will  be  an  indefinite  delay  in 
preparing  these  for  publication,  because  of  other  demands  on  his  time. 
Hence  the  project  has  been  terminated  for  the  time  being. 

Publication:  "Solution  of  the  telegrapher's  equation  with  boundary 

conditions  on  only  one  characteristic,"  by  G.  E.  Forsythe;  accepted  for 
publication  in  the  NBS  Journal  of  Research. 


THEORETICAL  AND  PRACTICAL  PROBLEMS  IN  THE  APPLICATION 
OF  A.D.C.M.  TO  THE  DETERMINATION  OF  ORBITS 
Project  11 .1/1-49-AM3 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  previous  status  report  for  this  project 

described  progress  on  a manuscript  entitled  "A  Modification  and  Ap- 
praisal of  Gibbs'  Orbit  Method,"  and  a trip  to  the  Yale  University  Libra- 
ry to  consult  some  recently  discovered  manuscripts  of  Willard  Gibbs  . 

As  a result  of  the  report  on  these  manuscripts  a grant  of  somewhat  over 
$1,000  was  made  by  the  National  Academy  of  Sciences  Watson  Fund  to  con- 
solidate the  new  material  in  the  manuscript  and  complete  the  work. 
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The  previous  report  also  described  progress  on  a text  on  rocket 
navigation.  In  the  present  quarter  two  new  chapters  were  completed;  one 
was  completely  revised;  and  several  were  completely  outlined.  Numerical 
investigations  and  examples  make  up  a sizable  portion  of  the  work  completed. 

Publications:  (l)  "Introduction  to  tables  for  rectilinear  and 

nearly-rectilinear  ('nearly  parabolic')  orbits,"  to  be  published  in  the 
NBS  Applied  Mathematics  Series;  see  also  project  11.1/2-48-3.  (2)  "A  mod- 

ification and  appraisal  of  Gibbs'  orbit  method,"  by  S.  Herrick;  IN  MANU- 
SCRIPT in  possession  of  the  author.  (3)  The  text  on  rocket  navigation 
cited  above . 


THEORETICAL  AMPLITUDE  DISTORTION  IN  A 
CONDENSER  MICROPHONE  CIRCUIT 
Project  11.1/1-49-AM4 

Origin:  Department  of  .Physics,  UCLA  Authorized  11/1/48 

Sponsor:  Office  of  Naval  Research,  USN  Terminated  9/30/49 

Manager:  G.  E.  Forsythe 

Objective:  ( l)  To  study  mathematically  the  relation  between  the 

variable  (input)  capacitance  of  the  condenser  element  and  the  variable 
(output)  voltage  across  the  resistance  of  an  RC  condenser  microphone  cir- 
cuit. In  particular,  to  find  conditions  under  which  the  voltage  faith- 
fully reproduces  the  capacitance,  as  the  resistance  becomes  infinite. 

(2)  To  calculate  the  amplitudes  of  the  various  harmonics  in 
the  output  voltage,  corresponding  to  a pure  sinusoidal  input  and  a finite 
resistance . 

Background:  A standard  reference  on  the  condenser  microphone, 

E.  C.  Wente,  Physical  Review  ID,  39-63  (1917)  makes  some  approximations 
which  raised  questions  about  the  practical  effect  of  dead  capacitance  in 
parallel  with  a condenser  microphone.  P.  M.  Morse,  "Vibration  and  sound," 
gives  a proof  which  is  valid  only  for  infinite  resistance.  It  is  be- 
lieved that  no  one  ha.3  presented  an  exact  formula  for  the  amplitudes  of 
the  various  harmonics  for  finite  R.  All  these  matters  are  dealt  with  in 
the  present  project  which  was  suggested  by  Professor  R.  W.  Leonard  of 
the  Department  of  Physics,  UCLA. 

Comments:  Professor  Leonard  plans  an  experiment  to  determine  the 

amplitudes  mentioned  under  (2)  above.  After  his  circuit  constants  are 
known,  a calculation  may  be  carried  out  to  confirm  his  results.  It  would 
involve  rapidly  convergent  series  of  Bessel  functions. 

Status:  TERMINATED.  The  project  has  been  inactive  for  some  time, 

pending  an  experiment  to  be  performed  by  Professor  Leonard.  Since  it 
now  appears  that  there  will  be  an  indefinite  delay  in  performing  this  ex- 
periment, the  project  has  been  terminated  for  the  present. 
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NUMERICAL  METHODS  IN  CONFORMAL  MAPPING 
Project  11.1/1-49-CM1 


Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  19^9  issue. 

Status:  CONTINUED.  This  project  was  perhaps  the  most  active  of  the 

research  projects  of  the  Institute  during  the  quarter.  The  activity  was 
largely  analytical  and  expository,  much  of  it  being  centered  about  the 
preparation  of  papers  for  the  "Proceedings"  of  the  Symposium  held  in  June 
(see  Publications).  However,  at  the  end  of  the  quarter  several  of  the 
project  managers  (notably  Drs . 'Warschawski  and  Ahlfors)  had  arrived  at 
points  in  their  work  at  which  they  felt  that  extensive  numerical  experi- 
mentation was  now  justifiable  and  desirable. 

Dr.  Beckenbach,  who  is  the  editor  of  the  "Proceedings,"  worked 
during  the  quarter  almost  exclusively  on  editorial  tasks  connected  with 
preparing  the  manuscripts  for  the  printer. 

Dr.  Seidel,  who  left  the  Institute  to  return  to  his  academic 
post  after  the  third  week  of  the  quarter,  wrote  up  his  Bibliography  for 
the  Proceedings  [[Publication  (lj.)J,  and  further  revised  his  translation 
Publication  (2)  . 

Dr.  Warschawski  spent  much  of  the  time  writing  up  his  paper  for 
the  Proceedings  [Publication  (l,ii)].  This  paper  has  now  become  a com- 
prehensive review  of  results  dealing  with  the  variation  of  the  mapping 
function  corresponding  to  deformations  in  the  mapped  region.  Both 
"nearly  circular"  and  arbitrary  regions  are  considered.  Emphasis  is  placed 
on  theorems  of  a quantitative  type  which  give  actual  estimates  for  the 
modulus  of  the  difference  between  the  mapping  functions  for  two  regions 
which  are  "close"  to  each  other.  A considerable  amount  of  foreground  re- 
search entered  into  the  preparation  of  the  paper. 

Dr.  Warschawski  also  began  a study  of  the  error  term  involved 
in  the  Bieberbach  method  of  constructing  mapping  functions  by  the  use  of 
the  Rayleigh -Ritz  principle.  He  expects  to  continue  this  work  at  the  Uni- 
versity of  Minnesota;  it  will  be  supplemented  by  numerical  experimentation. 

Dr.  Ostrowski  completed  two  papers  giving  the  results  of  his  study 
of  the  convergence  features  of  the  method  of  conformal  mapping  developed 
by  Theodorsen  and  Garrick.  The  contents  of  these  papers  had  been  presented 
by  Dr.  Ostrowski  at  the  Symposium  on  Conformal  Mapping  held  at  the  Institute 
for  Numerical  Analysis  in  June,  and  the  two  papers  will  appear  in  the 
Proceedings  [Publications  (l,iii)  and  (1,  iv)J.  The  results  contained  in 
the  previous  quarter's  status  report  are  included  in  the  first  paper  [Pub- 
lication (l,iii)J.  Also  included  in  this  paper  is  the  further  result  that 
the  convergence  of  the  sequence  of  the  derivatives  of  approximating  func- 
tions can  be  established  under  less  restrictive  conditions.  The  second 
paper  [Publication  (l,  iv)[j  deals  with  the  matrix  analogue  of  the  Theodorsen 
and  Garrick  method.  The  method  starts  from  the  approximating  trigonometric 
expression  for  the  singular  integral  involved  in  Theodorsen 's  and  Garrick's 
method,  and  gives  a number  of  properties  of  the  corresponding  matrices. 

These  matrices  are  singular  but  present  many  analogies  with  orthogonal 
matrices.  The  determination  of  their  characteristic  numbers  enables  one 
to  show  that  the  analogue  of  Warschawski ' s convergence  condition  cannot  be 
improved.  The  numerical  procedure  involved  is  appropriate  for  use  with 
high-speed  computing  machinery. 

Dr.  Weinstein,  with  the  help  of  Mr . R.  J.  Diamond,  one  of  the 
summer  fellows  at  the  Institute,  wrote  up  the  paper  which  he  had  presented 
at  the  Symposium  [Publication  (1,  v)]  during  a two-weeks'  visit  to  the  In- 
stitute . 

Dr.  Lars  Ahlfors  was  present  at  the  Institute  only  for  the  month 
of  August.  During  his  visit  he  prepared  the  report  referred  to  below  in 
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Publication  (3).  This  report  discusses  various  ways  of  constructing  the 
mapping  function  which  maps  a polygon  on  a circle  or  half  plane.  The 
starting  point  of  this  report  is  as  follows:  Let  the  polygon  lie  in  the 

z-plane  with  vertices  at  zk,  k = 1,2, ..  .n,  and  let  the  angles  be  denoted 
by  *k  7T . The  mapping  on  the  unit  circle  is  then  g^ven  by 

p'n  *k-l 

z = C J it  (w  - ak)  dw, 

J k=l 

where  ak  = eiw^  and  the  constant  C are  unknown  parameters.  Three  of  the 
uj  can  oe  chosen  arbitrarily,  and  C determines  only  the  size  of  the  polygon. 
However,  in  general  one  is  required  to  determine  the  mapping  function  which 
lets  w = 0 correspond  to  a prescribed  point  zQ  in  the  polygon.  It  is  there- 
fore better  to  make  this  normalization  and  treat  all  differences  - u>k 

as  unknown . Now 


*°k+l  " ^k  " 2 w uk(zo) 

where  uk  (the  harmonic  measure)  is  the  harmonic  function  ir  the  polygon 
which  is  equal  to  one  on  the  side  zk  zk+l  an^  zdro  on  the  rest  of  the 
boundary.  The  most  important  problem  is  therefore  to  compute  the  harmonic 
measures  of  the  sides.  If  this  is  done  for  a fixed  z0  the  inverse  function 
z(w)  will  be  determined  explicitly  (up  to  a constant  factor)  by  the 
Schwarz-Christoffel  formula.  On  the  other  hand,  if  two  harmonic  measures 
are  computed  for  variable  z it  is  easy  to  determine  the  direct  mapping  func- 
tion (w(z)  will  be  the  intersection  of  two  circular  arcs).  The  report 
therefore  concentrates  on  finding  a numerical  method  for  calculating  the 
harmonic  measure  of  a polygon  side. 

Two  quite  different  procedures  for  calculating  the  harmonic  meas- 
ure are  thus  proposed.  The  first  is  based  on  polynomials  possessing  ex- 
tremal properties  on  the  boundary  of  the  region  to  be  mapped  in  the  unit 
circle.  These  polynomials  are  closely  related  to  the  Szego  orthogonal 
polynomials.  An  advantage  of  this  method  is  that  it  yields  simple  explicit 
formulas,  but  the  computation  of  the  formulas  is  probably  prohibitively 
time-consuming,  even  for  a.d.c.m.  Dr.  Ahlfors  also  considers  briefly 
the  computing  method  proposed  by  Bergman  involving  areal  as  opposed  to  line 
integrals,  but  expresses  a preference  for  the  line  integrals.  The  second 
method  is  a generalization  of  Schwarz's  "alternating  method,"  or  method  of 
overlapping  regions  . The  first  step  is  to  cover  the  region  to  be  mapped 
by  a set  of  overlapping  subregions  for  which  the  Dirichlet  problem  is 
solvable.  It  Is  shown  that  the  problem  reduces  to  the  solution  of  an  in- 
tegral equation  of  the  following  type: 

u(x)  = J K(x,t)u(t)dt  + F(x), 

if 

where  F(x)  and  K(x,t)  are  explicitly  known  in  a given  case  (they  are  ob- 
tained from  the  assumption  that  the  Dirichlet  problem  is  solvable  for  each 
of  the  covering  subregions)  . The  report  then  discusses  iteration  methods 
of  solving  the  integral  equations.  Dr.  Ahlfors  believes  that  the  second 
method  outlined  above  is  more  promising  for  numerical  work,  and  recommends 
strongly  that  it  be  tested. 

Publications:  (l)  "Proceedings  of  a symposium  on  the  construction 

and  application  of  conformal  mapping"  edited  by  E.  F.  Beckenbach;  to  be 
published  by  NBS . The  volume  will  include  the  following  papers  written  in 
connection  with  this  project:  (i)  "A  bibliography  of  numerical  methods 

in  conformal  mapping,"  by  W.  Seidel.  (ii)  "On  conformal  mapping  of 
variable  regions,"  by  S.  E.  Warschawskl . (iii)  "On  the  convergence  of 
Theodorsen's  and  Garrick's  method  of  conformal  mapping,"  by  A.  M.  Ostrowskl. 
(iv)  "On  a discontinuous  analogue  of  Theodorsen's  and  Garriok's  method," 
by  A.  M.  Ostrowski . (v)  "On  the  Helmholtz  problem  of  conformal  represen- 

tation," by  Albert  Weinstein. 
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(2)  "Conformal  representation  of  simply -and  multiply-con- 
nected regions,"  by  L.  Kantorovitch  and  others;  translation  from  the 
Russian  by  W.  Seidel;  probably  to  be  published  in  NBS  Applied  Mathe- 
matics Series  .' 

(3)  "Numerical  methods  in  conformal  mapping,"  by  Lars  Ahlfors; 
IN  MANUSCRIPT,  awaiting  supplementary  numerical  work. 


MONOGRAPH  ON  NUMERICAL  SOLUTIONS  OF  DIFFERENTIAL  EQUATIONS 

Project  11.1/1-49-ODE1 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Manager:  W.  E.  Milne 

Objective:  To  prepare  a systematic  exposition  of  the  theory  and 

techniques  of  integrating  differential  equations,  with  special  emphasis 
on  ordinary  differential  equations  . The  treatise  will  include  chapters 
on  the  analysis  of  remainder  terms,  the  accumulation  of  errors,  relaxation 
methods,  and  other  topics  pertinent  to  a.d.c.m. 

Background:  The  literature  of  numerical  methods  in  differential 

equations  is  at  present  scattered,  incomplete,  and  often  aimed  at  very 
special  cases.  A convenient,  up-to-date  treatise  is  very  much  needed. 

(See  also  project  11.1/1-49-  DE2 . ) 

Comments:  Project  11 . 1/1-49-0DE2  is  intended  to  furnish  some  of  the 
groundwork  for  this  project. 

Status:  TERMINATED.  The  project  was  terminated  when  the  project 

manager  returned  to  his  regular  academic  post  at  Oregon  State  College. 

Work  is  continuing  there  on  the  project.  It  is  expected  that  Dr.  Milne 
will  return  to  the  staff  of  the  Institute  some  time  during  the  fiscal 
year  1951  j at  which  time  the  project  will  be  completed. 


Authorized  11/26/48 
Terminated  9/30/49 


STUDIES  IN  NUMERICAL  INTEGRATION  OF 
ORDINARY  DIFFERENTIAL  EQUATIONS 
Project  11.1/1-49-0DE2 


Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  1949  Issue. 

Status:  CONTINUED.  Dr.  Milne  left  the  I.N.A.  for  his  academic  post 

at  the  beginning  of  the  quarter,  but  a considerable  amount  of  work  con- 
tinued to  be  done  on  this  project,  both  at  the  I.N.A.  and  by  Dr.  Milne 
at  home . 

In  particular.  Publication  (3)  was  prepared  by  Dr.  Ostrowski 
and  Mr.  Jackson  from  correspondence  with  Professor  Bieberbach.  The  cor- 
respondence was  initiated  by  Dr.  Ostrowski  when  he  acted  as  a reader  for 
Publication  (2)  and  noticed  that  no  proof  was  available  for  a pertinent 
formula  in  Bieberbach 's  "Theorie  der  Dlff erentialgleichungen  . " Bieberbach' s 
reply  provided  an  Improvement  of  the  formula  but  contained  minor  algebraic 
errors,  which  were  corrected  by  Dr.  Ostrowski  and  Mr.  Jackson.  A princi- 
pal result  is  that  if  y(x)  denotes  the  solution  of 
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on  |x  - x | ^ a,  and  y]_(x)  denotes  the  approximation  to  y(x)  obtained  by  the 
Runge-Kutta  rule  (starting  with  correct  initial  values),  then 

|y(x)  - yi(x)|  K|x  xQ|5, 

where  K is  a constant  explicitly  given  in  terms  of  bounds  of  f(x,y)  and 
its  derivatives  so  chosen  as  to  be  governed  simultaneously  by  the  mag- 
nitude of  a. 

Dr.  Forsythe  studied  the  results  of  Drs . Huskey  and  Hartree 
(NBS  J.  Res.  ji2,  57  to  62  (1949)  on  numerical  integration  of  the  system 


dx 

dt 


J, 


x 


by  the  Heun  method.  Comparison  of  the  numerical  solutions  with  the  true 
sine-cosine  solution  of  the  system  revealed  surprising  deviations  of  the 
round-off  errors  from  those  expected.  In  the  article  cited,  Hartree 
gave  one  criterion  for  determining  certain  values  of  t where  qualitatively 
bad  round-off  errors  might  be  expected,  but  most  of  the  anomalies  were 
still  left  unexplained.  In  the  present  investigation  the  Hartree  criterion 
has  been  amplified  so  that  it  provides  an  estimate  of  the  round-off  error 
expected  to  occur  in  the  vicinity  of  each  point  t.  This  estimate  actually 
explains  most  of  the  anomalies  found  in  the  Huskey  data  and  gives  a method 
of  compensating  for  these  errors  to  a considerable  degree. 

Huskey's  round-off  errors  behaved  mysteriously  essentially 
because  the  mathematical  functions  involved  have  subtle  number-theoretical 
properties.  It  is  felt  that  too  much  ad  hoc  analysis  would  be  required  to 
compensate  every  Integration  for  its  round-off  errors.  Therefore,  in 
numerical  integration  high-speed  automatic  digital  computing  machines,  it 
is  proposed  that  one  force  the  round-off  to  obey  known  statistical  dis- 
tribution laws  by  the  following  device:  Whenever  a number  is  to  be  cur- 

tailed and  rounded,  let  the  rounding  be  "up"  or  "down"  according  to  some 
unbiased  random  procedure.  Experiments  on  IBM  machinery  are  testing  this 
method  on  the  above  system. 

Publications:  (1)  "A  note  on  the  numerical  integration  of  differential 

equations,"  by  W.  E.  Milne;  to  be  published  in  the  NBS  Journal  of  Research. 
(2)  "Note  on  the  Runge’-Kutta  method,"  by  W.  E.  Milne;  IN  MANUSCRIPT  in  pos- 
session of  the  author.  (3)  "On  the  remainder  of  the  Runge-Kutta  formula 
in  the  theory  of  ordinary  differential  equations,"  by  L.  Bieberbach; 
edited  and  revised  by  A.  M.  Ostrowski  and  L.  K.  Jackson  of  the  Institute 
for  Numerical  Analysis;  submitted  to  the  Quarterly  of  Applied  Mathematics. 

Note:  A paper  by  Dr.  Milne  entitled  "Numerical  determination 

of  characteristic  numbers,"  previously  listed  under  this  project,  now  ap- 
pears in  the  publications  for  project  11. 1/1-50-3- 
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APPROXIMATE  DETERMINATION  OF  THE  CHARACTERISTIC  FUNCTIONS  AND 
CHARACTERISTIC  VALUES  OF  ARBITRARY  LINEAR 
DIFFERENTIAL  OPERATORS 
Project  11.1/1-49-ODE4 

Origin:  NBS  Authorized  2/15/49 

Sponsor:  Office  of  Naval  Research,  USN  Terminated  9/30/49 

Manager:  C.  Lanczos 

Objective:  To  develop  a practical  and  economical  method  of  ob- 

taining the  characteristic  values  and  functions  of  arbitrary  linear 
differential  operators. 

Background:  In  many  problems  of  physics  and  engineering,  the 

determination  of  the  characteristic  solutions  of  a linear  differential 
equation  is  of  great  importance.  The  frequently  employed  Rayleigh-Ritz 
method  is  applicable  only  to  self-adjoint  systems,  it  requires  integra- 
tions, and  is  basically  restricted  to  the  characteristic  value  which 
is  least  in  absolute  value.  The  present  study  is  aimed  at  developing 
a method  which  will  not  have  these  limitations. 

Comments:  The  principle  underlying  the  present  investigation 

consists  in  constructing  a trial  function  which  together  with  its  iter- 
ations up  to  a certain  order  satisfies  the  given  boundary  conditions, 
and  then  determining  the  characteristic  equation  by  a method  entirely 
analogous  to  that  used  in  project  11 . l/l-49-AEl . 

Status:  TERMINATED.  It  was  decided  that  the  scope  of  the  objec- 
tive should  be  widened,  so  project  11.1/1-50-3  was  set  up  to  replace 
this  project.  See  current  status  report  for  project  11.1/1-50-3. 


STUDIES  IN  NUMERICAL  INTEGRATION  OF 
PARTIAL  DIFFERENTIAL  EQUATIONS 
Project  11.1/1-49-PD1 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Although  the  project  manager  returned  to  his 

academic  post,  some  numerical  work  continued  on  the  project  at  the 
Institute  and  the  project  manager  returned  to  the  staff  as  a consultant 
for  a short  period  in  September.  The  study  mentioned  in  the  previous 
status  report  was  completed  and  was  incorporated  into  a paper  presented 
at  the  Symposium  held  at  the  Harvard  Computation  Laboratory  in  Sep- 
tember . 

Publications:  "Numerical  methods  associated  with  Laplace's 

equation, "by  W.  E.  Milne;  IN  MANUSCRIPT  in  possession  of  the  author; 
to  appear  In  the  Proceedings  of  the  above  Symposium. 
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SAMPLING  TECHNIQUES  FOR  SOLVING  PARTIAL  DIFFERENTIAL 
EQUATIONS  AND  INTEGRO -DIFFERENTIAL  EQUATIONS 
Project  11.1/1-49-PM1 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Experimental  work  on  the  integration  of  Laplace's 

equation  was  performed  by  Mr.  Ernest  Elyash,  a summer  fellow  of  the  Insti- 
tute, under  the  direction  of  William  Feller.  The  joint  RAND-INA  Seminar 
on  Stochastic  Processes  continued  to  meet  regularly  every  Thursday  through- 
out the  quarter  under  the  direction  of  William  Feller,  who  spoke  at  a 
number  of  the  meetings.  Other  papers  were  given  by  Mr.  H.  Kahn  of  RAND, 
and  Mr.  Kac  of  the  Institute. 

Mr.  Feller  is  preparing  an  exposition  on  the  connection  between 
stochastic  processes  and  differential  equations,  based  on  his  lectures 
to  the  summer  fellows  (see  project  11.1/4-49-7)  and  to  the  joint  RAND-INA 
seminar . 

Publication:  "Proceedings  of  a Symposium  on  the  Monte  Carlo  Method," 

a volume  presenting  the  papers  given  at  the  Symposium  on  June  29,  30,  and 
July  1,  1949,  is  In  preparation. 


GENERATION  AND  TESTING  OF  RANDOM  DIGITS 
Project  11 .1/1-49-PM2 

Origin:  NBS 

Sponsor:  Office  of  Naval  Research,  USN 

Managers:  G.  E.  Forsythe,  R.  P.  Peterson,  Jr., 

E.  Yowell,  Jr. 

Objective:  (l)  To  devise  machine  methods  of  generating  random 

digits  by  arithmetical  processes.  (2)  To  test  whether  the  methods 
give  random  digits  suitable  for  use  in  various  Monte  Carlo  sampling 
processes . 

Background:  Random  digits  (or  random -sampling  digits)  are  needed 

for  many  purposes  in  "experimental  mathematics,"  applied  statistics,  and 
numerical  analysis.  (See  for  example  project  11 .1/1-49-PM1 . ) Tippett 
and  Kendall-Smith  in  England  have  published  tables  of  random -sampling 
digits  obtained  by  laborious  hand  work  with  tables  and  machines  of  the 
roulette-wheel  type.  The  RAND  Corporation  has  generated  random  digits 
from  the  thermionic  noise.  It  is  inherent  in  the  nature  of  a.d.c.m., 
however,  that  neither  precomputed  digits  now  produced  by  a non-reprodu- 
cible process  are  very  well  suited  for  calculations.  It  therefore  seems 
essential  to  have  techniques  for  internal  generation  of  random  digits 
(l)  which  will  be  repeatable,  and  (2)  which  will  give  digits  which  be- 
have like  the  products  of  a truly  random  process. 

Comments:  Some  of  the  atomic  Energy  Commission  Laboratories  have 

experimented  with  the  arithmetic  generation  of  sampling  digits.  However, 
the  reports  at  the  Symposium  on  the  Monte  Carlo  Method  held  at  the  INA 
on  June  29,  30,  July  1,  1949,  showed  a need  for  more  basic  Investigations. 


Authorized  4/1/49 
Completed  9/30/49 
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Status:  COMPLETED. 

Publication:  The  results  will  appear  In  the  "Proceedings  of  a 

Symposium  on  the  Monte  Carlo  Method,"  to  be  published  by  the  NBS . 


DETERMINANTS  OF  ORTHOGONAL  POLYNOMIALS 
Project  11.1/1-49-SF1 

Origin:  NBS  Authorized  11/29/48 
Sponsor:  Office  of  Naval  Research,  USN  Completed  9/30/49 
Managers:  E.  F.  Beckenbach,  W.  Seidel,  G.E. Forsythe 


Objective:  To  derive  identities,  inequalities,  and  convexity 

properties  of  minors  of  recurrent  determinants  whose  elements  are 
orthogonal  polynomials,  e.g.. 


A (n>  h<  k;  x)  = 


Pn<x> 


Pn  , (x) 

n+h  ' 


n+k 


(x)  P . , (x) 
v ' n+h+k  ' 


where  Pn(x)  is  the  Legendre  polynomial  of  degree  n • and  to  develop  ap- 
plications of  these  results/  especially  to  the  problem  of  moments. 


Background:  This  Is  primarily  a long-range  project  aimed  at 

extending  general  mathematical  knowledge.  However,  orthogonal  poly- 
nomials are  of  great  importance  in  applied  mathematics,  and  any 
further  information  concerning  them  has  an  appreciable  probability 
of  being  immediately  useful . 


Status:  COMPLETED.  In  a forthcoming  paper  [Publication  (2)J, 

E.  F.  Beckenbach,  W.  Seidel  and  0.  Szasz  have  studied  recurrent  de- 
terminants of  ultrasphericai  polynomials  and  have  shown,  among  other 
things,  that  the  quadratic  forms 


(A) 


U 

V 


Z PU)(x)/PU)(l)  u 

0 [_  /*.YK  J /*■ 

where  Pn''*/(x)  is  the  ultrasphericai  polynomial  of  degree 
are  positive  definite  for  A > o,  n = 0,1,2,  ... , and  x > 1 . 
present  paper  it  is  shown  that  the  Hermitian  forms 


n and  order  A , 
In  the 


I 

L 

are  positive  definite  for  A A 0,  n = 0,1,2,...  and  -1<x<1,  a result 
conjectured  by  G.  Szego.  (This  is  the  problem  discussed  in  the  second 


paragraph  of  the  previous  status  report  for  this  project 
suit  follows  from  the  fact  that  the  harmonic  function. 


This  re- 


jr  + R 


(f  [pn(A>Cx)/Pn(A>  (1)]  zn] 


is  positive  in  the  unit  circle  Jzl<l  for  -l^x<l  and  A > 0 . For  the 
special  case  A = 1 this  result  was  proved  earlier  by  G.  Szego. 


Publications:  (l)  "Second  order  determinants  of  Legendre  poly- 
nomials," by  G.  E.  Forsythe;  to  be  submitted  to  the  Duke  Mathematical 
Journal.  (2)  "Recurrent  determinants  of  orthogonal  polynomials,  Part 
I:  Legendre  and  ultrasphericai  polynomials,"  by  E.  F.  Beckenbach, 
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W.  Seidel,  and  0.  Szasz;  accepted  by  the  Duke  Mathematical  Journal, 
subject  to  revision  and  condensation. 


DETERMINATION  OP  EXTREMALS  OF  FUNCTIONALS 
Project  11. 1/1-50-1 


Origin:  NBS  Authorized  7/1/49 

Sponsor:  Office  of  Naval  Research,  USN 

Manager:  M.  R.  Hestenes 

Objective:  To  study  numerical  and  analytical  methods  of  deter- 

mining the  maxima  and  minima  of  functionals;  in  particular,  integrals 
of  the  types  customarily  studied  in  the  calculus  of  variation. 

Background:  A substantial  part  of  modern  physical  theory  can  be 

stated  in  terms  of  variational  principles.  In  view  of  this  fact,  the 
calculus  of  variation  is  playing  a continually  expanding  role  in  the  sol- 
ution of  problems  in  applied  physics.  The  existing  literature  on  cal- 
culus of  variation  Is  extensive,  but  the  numerical  aspects  have  not  kept 
pace  with  the  theoretical  ones. 

Status:  NEW.  Publication  (l)  was  completed,  and  a certain  amount 

of  time  was  spent  in  revising  Publication  (2). 

Publications:  (l)  "Some  elementary  problems  in  the  calculus  of 

variations,"  by  M.  R.  Hestenes;  submitted  to  the  Mathematics  Magazine, 
also  to  be  published  as  a chapter  in  a survey  of  mathematics  edited 
by  Professor  Glenn  James  of  U.C.L.A.  (2)  "Quadratic  forms  In  Hilbert 
space,  with  applications  in  the  calculus  of  variation,"  by  M.  R.  Hestenes; 
accepted,  subject  to  revision,  by  the  American  Journal  of  Mathematics. 


DETERMINATION  OF  CRITICAL  POINTS 
Project  11.1/1-50-2 

Origin:  NBS  Authorized  7/1/49 

Sponsor:  Office  of  Naval  Research,  USN 

Manager:  M.  R.  Hestenes 

Objective:  To  study  numerical  methods  of  determining  critical 

points  of  functions  of  more  than  one  variable,  and  in  particular,  to 
develop  practical  methods  for  locating  maxima  and  minima  of  functions 
of  many  variables . 

Background:  Critical  point  determinations  arise  in  the  theory  of 

aircraft  design,  the  theory  of  games,  and  in  many  other  physical  and 
economic  applications . 

Status:  NEW.  Four  closely  related  methods  of  locating  maxima 

and  minima  of  functions  were  studied.  Each  one  is  a generalization 
or  modification  of  the  Newton -Raphson  iterative  method  of  finding  roots 
of  equations.  In  the  first  method,  first  and  second  order  derivatives 
are  used  to  find  the  roots  of  the  first  derivatives,  exactly  as  in 
the  Newton -Raphson  method.  The  second  method  is  really  only  a geo- 
metric interpretation  of  the  first  and  third  methods.  In  this  method 
level  surfaces  of  the  function  in  the  neighborhood  of  a critical  point 
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are  approximated  by  ellipsoids.  The  centers  of  these  ellipsoids  are 
found,  and  a process  is  set  up  whereby  the  successive  centers  converge 
toward  the  critical  point  in  question.  In  the  third  method,  the  function 
is  interpolated  by  successive  polynomials  of  the  second  degree,  and  the 
respective  minifna  or  maxima  of  these  polynomials  are  then  determined. 

The  last  method  is  the  so-called  method  of  steepest  ascent.  In  this 
method  the  first  problem  is  that  of  selecting  a norm  whereby  to  give  a 
meaning  to  the  term  "steepest  ascent."  It  is  found  that  when  a suitable 
norm  has  been  chosen,  the  method  reduces  simply  to  the  first  one  men- 
tioned above.  The  first  three  methods  are  about  equally  satisfactory 
as  regards  convergence;  choice  will  depend  on  convenience  in  specific 
cases.  Extensions  to  the  calculus  of  variations  have  been  considered. 
Progress  on  this  aspect  will  be  reported  on  in  later  status  reports  for 
project  ll.l/l-50-l.  Random  sampling  methods  have  been  investigated 
in  a preliminary  way,  but  without  satisfactory  results  as  yet. 


CALCULATION  OF  EIGENVALUES,  EIGENVECTORS,  AND  EIGENFUNCTIONS 

OF  LINEAR  OPERATORS 
Project  11.1/1-50-3 

Origin:  NBS  Authorized  7/1/49 

Sponsor:  Office  of  Naval  Research,  USN 

Managers:  C.  Lanczos,  W.  Feller,  G.  E.  Forsythe,  M.  Kac 

Objective:  To  develop  practical  and  economical  methods  of  solving 

eigenvalue  problems,  to  be  used  in  connection  with  various  types  of 
computing  equipment . 

Background:  The  determination  of  the  eigenvalues  (or  characteristic 

values),  eigepvectors , and  eigenfunction  expansions  of  matrices,  dif- 
ferential equations,  integral  equations,  etc.  is  of  fundamental  impor- 
tance in  many  problems  of  physics  and  engineering.  The  interest  of  the 
Institute  for  Numerical  Analysis  in  such  problems  originated  in  the  fact 
that  at  the  time  when  work  first  began  on  the  predecessor  of  this 
project  (see  Comments  below),  current  methods  of  determining  the  eigen- 
values of  matrices  were  inadequate  for  certain  matrices  which  arise  in 
vibration  analysis  problems . 

Comments:  This  project  represents  a continuation  and  amplifica- 

tion of  projects  11 .1/1-49-AE1  and  11 . 1/1-49-ODE4 . 

Status:  NEW.  (I)  The  manuscript  of  Dr.  Lanczos  referred  to  in 

the  previous  quarter's  status  report  for  project  11 . 1/1-49-AE1  was 
incorporated  into  the  more  general  treatment  represented  below  in  Pub- 
lication (l).  This  treatment  also  includes  the  material  on  integral 
equations  and  differential  operators  described  in  earlier  status  re- 
ports for  project  11 . 1/1-49-0DE4 . As  certain  mathematical  details  of 
the  Lanczos  method  have  already  been  given  in  these  earlier  status  re- 
ports, no  attempt  will  be  made  here  to  describe  the  method  mathematically 
as  it  finally  emerges  in  Publication  (l).  The  advantages  claimed  for 
the  method  are  as  follows:  (a)  The  iterations  are  used  in  the  most 

economical  fashion,  obtaining  an  arbitrary  number  of  eigenvalues  and 
eigensolutions  by  one  single  ,set  of  iterations,  without  reducing  the 
order  of  the  matrix.  (b)  The  rapid  accumulation  of  fatal  rounding  errors, 
common  to  all  iteration  processes  if  applied  to  matrices  of  high  dispersion 
(large  "spread"  of  the  eigenvalues),  is  effectively  counteracted  by  the 
method  of  "minimized  iterations."  (c)  The  method  is  directly  translatable 
into  analytical  terms  by  replacing  summation  by  integration.  We  then  get 
a rapidly  convergent  analytical  iteration  process  by  which  the  eigenvalues 
and  eigensolutions  of  linear  differential  and  integral  equations  may  be 
obtained . 

(II)  In  the  course  of  studying  and  testing  the  Lanczos  method 
referred  to  above.  Dr.  Milne  developed  a procedure  for  calculating  eigen- 
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values,  based  essentially  on  the  Lanczos  method,  but  modified  in  such  a 
way  as  to  provide  a simple  numerical  routine  for  the  computation . In 
Publication  (2),  this  procedure  is  set  forth, < together  with  a number  of 
numerical  illustrations  . The  exposition  is  limited  to  second  order  dif- 
ferential equations.  The  essential  steps  in  Dr.  Milne's  procedure  are 
these : 

z 

Let  L(u)  = £ Pv(x)  D2“k  u 
C K 

be  the  differential  operator  in  question.  The  problem  is  to  find  the 
characteristic  values  of  for  which  the  differential  system 

L(u)  + A 2u  = 0 , 

u + g Du  = 0 at  x = a , 

u + G Du  = 0 at  x = b , 

possesses  non-zero  solutions  in  a^x^b.  A fictitious  partial  differ- 
ential system  Is  now  introduced; 

V = L(V) 

V + g Dx  V = 0 at  x = a 

V + G Dx  V = 0 at  x = b 

V = F(x)  when  t = 0 
Dt  V = 0 when  t = 0 

A solution  in  the  form 

(*)  V = Y-  ckuk(x)  cosAkt 

is  assumed.  A network  is  set  up  in  the  (x,t)  space  and  the  partial  dif- 
ferential system  Is  replaced  by  an  approximately  equivalent  difference 
equation.  The  values  of  the  function  V at  the  vertices  of  the  network 
may  be  considered  as  coefficients  of  a set  of  homogeneous  simultaneous 
linear  equations  with  a non-trivial  solution,  and  by  substitution  of  (*) 
into  the  equation  a corresponding  set  of  equations  involving  the  eigenval- 
ues Ak  is  derived.  From  these  It  is  evident  that  the  eigenvalues  of  the 

difference  equations  are  proportional  to  the  roots  of  an  equation  of  the 

type  n 

21  Ak  cos  k fL>-  = 0 . 

(Ill)  Dr.  W.  Feller  and  Dr.  G.  E.  Forsythe  collaborated  on  a 
study  of  the  transformation  of  a given  square  matrix  A Into  another  matrix 
B such  that  the  characteristic  values  and  vectors  of  A can  be  determined 
from  those  of  B.  The  initial  ideas  were  mainly  those  of  Dr.  Feller.  As 
of  the  end  of  the  quarter,  a joint  paper  was  being  prepared  [Publication 
(5)  below].  Technical  details  of  this  work  are  as  follows: 

Two  well-known  but  previously  unrelated  numerical  procedures 
for  obtaining  the  eigenvalues  and  eigenvectors  of  an  n-th  order  square 
matrix  A both  begin  with  an  iteration  to  get  one  characteristic  value  A ]_ 
and  its  associated  vectors,  followed  by  a transformation  of  A to  a new 
matrix  B from  which  the  characteristic  value  A^  is  absent.  After  ob- 
taining by  iteration  one  characteristic  value  of  B and  associated  vectors, 
one  can  reconstruct  another  characteristic  value  and  vectors  of  A.  One 
can  then  start  over  with  B.  (Lanczos'  procedure  referred  to  in  (I)  and 
(II)  above  is  of  another  type.)  One_  procedure,  "deflation,"  is  symmetry- 
preserving (meaning  that  if  A Is  symmetric,  then  so  is  B)  , but  not  order- 
reducing  (order-reducing  means  that  B is  of  order  n-l) . The  other  proce- 
dure, the  Frazer-Duncan -Collar  reduction  method,  is  order-reducing,  but 
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not  symmetry-preserving. 

In  the  present  investigvation  a two  dimensional  set  S of 
matrix  transformations  of  A to  B is  found,  each  member  of  which  permits 
one  to  proceed  from  B's  characteristic  values  and  vectors  to  these  of  A. 

The  deflation  procedure  and  the  Frazer-Duncan -Collar  reduction  procedure 
both  appear  as  natural  special  members  of  the  set  S.  A certain  one- 
dimensional subset  of  S is  symmetry -preserving,  while  another  one-dimen- 
sional subset  is  order-reducing.  Two  special  transformations  are  both 
order ■ reducing  and  symmetry-preserving  and  are  therefore  recommended  for 
practical  work. 

The  theorems  apply  to  both  defective  and  non-defective  matrices. 

(IV)  Dr.  M.  D.  Donsker  and  Dr.  M.  Kac  developed  a method  for 
finding  the  lowest  eigenvalue  and  corresponding  eigenfunctions  of  Schro- 
dinger ' s equation  by  the  so-called  Monte  Carlo  technique.  Their  work  is 
contained  in  the  paper  cited  in  Publication  (4)  . Their  method  is  based 
on  the  following  considerations: 

Let  X,,  X2,  X^ , . . . be  independent  identically  distributed  ran- 
dom variables  each  having  mean  0 and  standard  deviation  1 and  let 
Sj.  - X^  + X2  + . . .+  X^ . Under  certain  general  assumptions  on  V(x),  the 
most  severe  of  which  is  that  V(x)  be  non -negative , it  can  be  shown  that 
the  limiting  distribution  function  <r  (cx,t)  of  the  random  variable 


n k<nt  JxT 


is  such  that 
(*) 

where  ^ (x)  is 


oo  OO  OO 

J J e~c<~3t  d^  cr  (ot,  t)dt  = / ^ (x)dx 

0 0 -OO 

the  fundamental  solution  of  the  differential  equation 


1 

2 


subject  to  the  conditions 


( s + V(x) ) p = 0 


f (x)  0 


x ± oo 


j ^ ' ( x ) | < M x'^0 

J^'(+  0)  - {<,'(-0)  = 2. 


The  fundamental  solution  ^ (x)  of  the  aforementioned  differen- 
tial equation  is  expressible  in  terms  of  the  normalized  eigenfunction 
{^/.•(x)}  and  eigenvalues  Aj  of  the  one  dimensional  Schrodinger  eigen- 
value problem 

1 di£_  _ v(x)  fj  (x)  = _ x f,  A>0 

^ dx^ 


a s 
(**) 


</j  (x) 


Substitution  of  (**)  into  (*)  and  subsequent  inversion  with 
respect  to  s permits  the  lowest  eigenvalue  to  be  expressed  as 


A = lim  — i log  f e-<*  d cr(«x,t)  . 
1 £ -poo  t \J  ** 
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The  distribution  function  <r(<*,t)  is  then  calculated  by  a sampling  process 
in  which  the  value  of  S above  is  repeatedly  determined  for  a suitably  large 
value  of  n.  Generalizations  to  multidimensional  Schrodinger  equations  are 
considered.  The  theory  for  the  second  eigenvalue  is  also  stated. 

In  a numerical  example  in  which  Ap  was  actually  0.71,  with 
n = 200  the  Monte  Carlo  approximation  was  0.f4. 

(V)  Dr.  Kac  considered  the  following  classical  eigenvalue  problem: 
Let  B be  a simply  connected  two-dimensional  region  of  area  |Bl  bounded  by 
the  curve  C.  Within  B let  u = u(x,y)  satisfy  the  equation 

Au  + Aa  = 0,  with  u = 0 on  C. 

Then  according  to  a well-known  result  of  Weyl,  the  number  of  eigenvalues 
less  than  a given  A is  asymptotically  equal  to  / B I A /4t t as  A 00  . There 
is  a corresponding  result  due  to  Carleman  for  the  squares  of  the  eigen- 
functions . Dr.  Kac  found  new  short  proofs  of  these  results  based  on  the 
theory  of  diffusion. 

Publications:  (l)  "An  iteration  method  for  the  solution  of  the  eigen- 

value problem  of  linear  differential  and  integral  operators,"  by  C.Lanczos; 
to  be  submitted  to  the  NBS  Journal  of  Research.  (2)  "Numerical  determina- 
tion of  characteristic  numbers,"  by  W.  E.  Milne ; to  be  published  in  the 
NBS  Journal  of  Research.-  (3)  A paper  embodying  results  of  W.  Feller  and 
G.  E.  Forsythe  to  be  submitted  to  the  Quarterly  of  Applied  Mathematics. 

(4)  "A  sampling  method  for  determining  the  lowest  eigenvalue  and  the  prin- 
cipal eigenfunction  of  Schrodinger ' s equation,"  by  M.  D.  Donsker  and 
M.  Kac;  to  be  published  in  the  NBS  Journal  of  Research. 


MISCELLANEOUS  STUDIES  IN  PURE  MATHEMATICS 
Project  11.1/1-50-4 


Origin:  NBS  Authorized  7/1/49 

Sponsor:  Office  of  Naval  Research,  USN 

Managers:  Various  members  of  the  staff 

Objective:  To  conduct  various  studies  in  those  fields  of  pure 

mathematics  on  which  the  theory  and  practice  of  numerical  analysis  are 
based;  also  to  carry  out  certain  investigations  in  pure  mathematics  which  re- 
quire large  amounts  of  computing  and  so  may  result  incidentally  in  a 
contribution  to  the  science  of  computing. 

Background:  The  principal  scientific  purpose  of  this  project  is 

simply  the  general  expansion  of  mathematical  knowledge.  It  is  important 
in  conducting  research  work  at  the  scientific  level  of  that  of  the 
Institute  to  insure  freedom  of  inquiry  to  the  staff,  and  in  particular, 
to  avoid  setting  up  barriers  to  research  at  the  apparent  (and  probably 
quite  ephemeral)  boundaries  of  a specific  program.  However,  there  are 
certain  secondary  aims  of  this  project  which  make  it  contribute  tangibly 
to  the  program  of  the  Institute.  Two  of  these  may  be  inferred  directly 
from  the  statement  of  the  Objective.  A third  is  that  it  is  to  the  advan- 
tage of  the  Institute  to  have  mathematical  members  of  its  research  staff 
retain  their  proficiency  in  the  techniques  and  logic  of  pure  mathematics, 
as  effective  research  in  numerical  analysis  and  applied  mathematics  nearly 
always  makes  use  of  the  methods  of  pure  mathematics. 

Status:  NEW.  The  work  on  this  project  during  the  quarter  was 

done  mainly  by  Dr.  Otto  Szasz . He  wrote  two  papers,  Publications  (l)  and 
(2)  below,  on  the  Gibbs  phenomenon  as  it  appears  in  applying  certain  sum- 
mation methods  to  trigonometric  series,  and  a third  one,  Publication  (3), 
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on  the  Perron  method,  of  summability . 

A principal  result  of  the  first  paper  is  that  if  sn(t)  denotes 
the  Euler  mean  of  the  series  2 sin  nt  /n,  with  parameter  r lying  in  the 
interval 

r x 


0<r4  1,  then  lim  sn(tn)  =J  — y % dy 

o ® 

as  nt  -e  t , nt^-^0. 

The  second  paper  contains  a discussion  of  the  Gibbs  phenomenon 
for  Hausdorff  means.  Hausdorff  means  of  a sequence  {sn}  are  defined  by 

hn  = L (2)  sv  J u(l-u)n_v  d}A(u), 

where  f (u)  is  of  bounded  variation,  ^ (l)  =1,  j>  (+0)  = ^ (o)  =0. 

The  main  problem  here  is  to  give  conditions  on  the  function 
^ (u)  which  are  necessary  and  sufficient  in  order  that  the  hn  means  of 
the  series 


£ 


sin  nt 
n 


.ir  - t)  present  a Gibbs  phenomenon.  When  applied  to  Holder 


means  of  order  p > 0,  in  which  case 

f <u)  ’ HpT  / (log  ax’ 


the  existence  of  a constant  p0  is  proved, so  that  the  Gibbs  phenomenon 
occurs  for  p<pQ,  but  not  for  p > p0 . It  is  found  that  p = 0.5826...  . 

In  contrast  to  the  equivalent  Cesll.ro  means  Gronwall  found  the  critical 
value  0 . 4395 . . . . 

The  third  paper  contains  a comparison  of  the  Perron  and  Cesaro 
methods  of  summability.  The  Perron  transform  of  a series  Z an  is 

on  ^ 


^(x) 


Z a x x+1 
I n X+X  X+X+l 


x+n-1 

x+^+n-l 


The  series  is  summable  to  s in  the  Perron  sense  if  lim  0^(x)  = 3- 
main  result  of  this  paper  is  that  this  summability  method  is  more  powerful 
than  Ces&ro  summability  of  order  X-  The  essential  tool  is  an  integral 
representation  of  0^(x). 

Publications:  (l)  "On  the  Gibbs  phenomenon  for  Euler  means,"  by 

0.  Szasz;  accepted  for  publication  in  Acta  Litterarum  ac . Sclent.,  Szeged. 

(2)  "The  Gibbs  phenomenon  for  Hausdorff  means,"  by  0.  Szasz;  IN  MANUSCRIPT. 

(3)  "On  a summation  method  of  0.  Perron,"  by  0.  Szasz;  to  be  submitted  for 
publication  in  Mathematische  Zeitschrift. 
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2 . Mathematical  Tables 


MATHIEU  FUNCTIONS  II 
Project  11.1/  2-45-1 

Origin:  Applied  Mathematics  Panel,  NDRC 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  coefficients  De^  and  Dop  which  were  trans- 

mitted from  the  Computation  Laboratory  have  been  checked  in  connection 
with  a problem  for  the  Air  Materiel  Command.  The  computations  of  the 
functions  which  had  been  started  will  be  continued. 

Publication:  To  be  submitted  for  publication  in  the  NBS  Applied 

Mathematics  Series. 


DECOMPOSITIONS  OF  ARC  TANGENTS 
Project  11.1/2-48-1 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full ’pro ject  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  manuscript  is  being  typed,  and  the  final 

proofs  are  being  subjected  to  tests  under  the  supervision  of  J.  Todd. 

Publications:  The  table  will  be  submitted  for  publication  in  the 

NBS  Applied  Mathematics  Series;  the  table  will  Include  the  decompositions 
for  integral  arguments  n 4 2,100.  "A  problem  on  arc  tangent  relations," 
by  J.  Todd  has  been  accepted  by  the  American  Mathematical  Monthly. 


SPECIAL  TABLE  OF  BESSEL  FUNCTIONS 
Project  11.1/2-48-2 

Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Jan -Mar  1949  issue. 

Status:  CONTINUED.  Although  this  project  was  reported  completed, 

work  is  continuing  on  the  preparation  of  the  manuscript  for  publication. 

Publication:  To  be  submitted  for  publication  in  the  NBS  Applied 

Mathematics  Series. 
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ROCKET  NAVIGATION  TABLES 
Project  11.1/2-48-3 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC , USAF 

Full  project  description  appears  In  Jan-Mar  1949  Issue. 

Status:  CONTINUED.  Although  this  project  was  reported  completed, 

work  Is  continuing  on  the  preparation  of  the  manuscript  for  publication. 

Publication:  To  be  submitted  for  publication  in  the  NBS  Applied 

Mathematics  Series. 


TABLES  OF  E-j/z),  SECOND  QUADRANT 
Project  11.1/2-49-1 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue  . 

Status:  CONTINUED.  The  auxiliary  functions  Un(z)  and  Tn(z)  were 

completed  for  n4  17. 


PUNCHED  CARD  LIBRARY 
Project  11.1/2-49-2 

Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Comments:  A catalog  of  tables  on  punched  cards  which  are  on  file  at 

the  Institute  may  be  obtained  by  addressing  the  Institute  for  Numerical 
Analysis,  405  Hilgard  Avenue,  Los  Angeles,  24,  California.  Within  the 
limits  of  the  program  of  the  computation  unit  of  the  Institute,  tables 
will  be  duplicated  upon  request,  provided  the  requester  furnishes  the 
blank  cards.  Requests  should  be  addressed  directly  to  the  Institute. 

Status:  CONTINUED.  The  following  tables  on  punched  cards  were 

added  to  the  library: 

1.  sin  x,  cos  x , A , <5*  ; x = [o(. 0001)1;  11d]  with  maximum  erpor 
of  1 unit  in  last  place.  These  were  computed  in  connection 
with  the  study  of  round-off  error  by  G.  E.  Forsythe. 

2.  Table  of  Bessel-Clifford  Functions,  The  functions  _ 

_ n _q  _ n 

x ^ Jn(  2 v/x)  , x 2 In(2  v/x)  , x 2 Kn(2  v/x)  > x 2 In^x) 

are  now  available  on  punched  cards  for  x ranging  between  0 
and  410  at  various  Intervals,  for  n = 0 and  1.  For  x>6.2, 
the  functions 

2 v/x  -in  -2  v/x  -in 

e x K (2  v/x)  and  e x 1(2  A)  are  given. 
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3-  E2(x)  = Jt  (e'^/u^du,  ,x  = [0(  .001)2;  7Dj 

x = [2(  . 001)10;  5S  or  6s]  , x = [l0(  .001)16;  13D]  . 

4.  The  NBS  "Table  of  Spherical  Bessel  Functions,"  Vol. 

1,  has  been  punched,  but  not  checked.  North  American 
Aviation  Corp.  is  punching  the  NBS  "Table  of  Cir- 
cular and  Hyperbolic  Tangents  and  Cotangents."  The 
Consolidated  Vultee  Aircraft  Corp.  has  undertaken  the 
punching  of  the  NBS  "Tables  of  Natural  Logarithms," 
Vol.  Ill  and  IV. 


3 . Development  of  Automatic  Computing  Machinery 


AIR  MATERIEL  COMMAND  COMPUTING  MACHINE 
Project  11.1/22-49-1 

Origin:  Air  Materiel  Command,  USAF 

Sponsor:  Air  Materiel  Command,  USAF 

Full  project  description  appears  In  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  arithmetic  unit  contract  has  been  given  to 

the  Cole  Instrument  Company  of  Los  Angeles.  They  have  produced  prototypes 
of  both  the  M and  A Registers,  and  have  returned  them  to  be  checked.  The 
M Register  was  checked  and  found  satisfactory;  production  on  the  M Register 
is  now  under  way  at  the  Cole  Instrument  Company.  Checking  on  the  A Regis- 
ter will  soon  be  completed.  It  will  then  be  returned  to  the  Cole  Instru- 
ment Company  and  production  on  it  will  begin  at  that  time.  Delivery  of 
the  completed  chassis  for  the  arithmetic  unit  is  expected  by  not  later 
than  November  15 . 

The  control  unit  chassis  are  nearing  completion,  and  the  racks 
to  hold  these  chassis  are  in  the  process  of  being  manufactured.  Comple- 
tion of  the  racks  and  partial  checking  out  of  the  control  unit  are  anti- 
cipated by  November  1. 

A complete  cathode  ray  tube  memory  unit  has  been  built  and  oper- 
ated successfully  in  storing  512  points  on  the  face  of  the  tube.  Although 

this  unit  operated  entirely  satisfactorily  it  was  found  that  the  gating 
circuitry  could  be  simplified.  Hence,  a second  model  was  constructed 
making  use  of  simplified  gating  circuitry.  This  second  memory  tube  unit 
is  now  being  tested  in  the  cathode  ray  tube  test  unit. 

The  machine  will  be  put  into  operation  initially  with  only  256 
spots  on  the  face  of  the  cathode  ray  tubes,  - that  is,  256  words  in  the 
high  speed  memory.  However,  tests  with  the  present  memory  units  indicate 

that  there  would  be  no  great  difficulty  in  storing  1024  words  in  the  mem- 

ory with,  possibly,  a rule  about  reference  to  nearby  points  on  the  memory. 
(Experiments  by  Division  13  of  NBS  and  by  other  groups  have  indicated  that 
such  a rule  can  not  be  avoided)  . Coding  personnel  definitely  prefer  the 
larger  memory  even  at  the  price  of  such  a rule.  The  Zephyr  is  being  con- 
structed in  such  a way  that  the  memory  tube  capacity  can  be  increased  to 
1024  at  any  future  date  without  greatly  altering  the  rest  of  the  machine. 
From  the  above  it  will  be  seen  that  some  sacrifices  are  being  made  in  order 
to  complete  the  machine  as  quickly  as  possible. 
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LOGICAL  NOTATION  AND  BLOCK  DIAGRAM  SYMBOLISM  FOR  A.D.C.M. 

Project  11.1/22-49-2 


Origin:  NBS 

Sponsor:  Air  Materiel  Command,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue  . 

Status:  CONTINUED.  Comments  on  the  existing  terminology  and 

notation  have  been  received  from  various  interested  persons  in  the  field 
of  automatic  computing  machines  . When  deemed  worthwhile  by  the  committee 
these  comments  and  suggestions  have  been  incorporated  in  the  practices 
of  the  coding  group  at  the  Institute,  and  in  the  reports  issued  on  the 
automatic  computer  being  designed  and  constructed  under  project  11.1/22-49-1. 
It  is  hoped  that  in  this  manner  a system  of  terminology,  notation,  and 
symbolism,  which  will  be  widely  acceptable,  can  be  established.  The  sys- 
tem is  still  expanding,  and  is  continuously  being  subjected  to  practical 
test  by  the  coding  group.  It  Is  expected  that  the  period  of  test  will 
extend  into  the  early  stages  of  operation  of  the  machine  being  constructed 
in  project  11.1/22-49-1. 


SEMI-AUTOMATIC  INSTRUCTION  FOR  ELECTRONIC  DIGITAL  COMPUTERS 
Project  11.1/22-50-1 


Origin:  NBS  Authorized  7/1/49 

Sponsor:  Air  Materiel  Command,  USAF 

Manager:  H.  D.  Huskey 

Objective:  To  establish  a system  of  semi-automatic  instruction  for 

all  general  purpose  electronic  digital  computing  machines  using  as  a 
model  the  automatic  computing  machine  being  constructed  under  project 
11.2/22-49-1. 

Background:  With  the  advent  of  general  purpose  high-speed  digital 

computing  machines  the  need  for  a simplified  system  of  coding  has  become 
evident.  Such  a system,  by  making  the  machine  do  much  of  the  tedious  and 
repetitious  part  of  coding,  would  greatly  decrease  the  work  of  the  coding 
staff . 


Status:  NEW.  Dr.  Huskey  presented  his  preliminary  results  in  a 

talk  entitled  "Semi-automatic  instruction  on  the  Zephyr"  at  a Symposium 
on  Large-Scale  Calculating  Machinery  which  was  held  at  Harvard  University 
in  September  1949-  The  subject  has  also  been  discussed  with  Mr.  E.F. Moore, 
who  has  submitted  a paper  on  the  subject  to  the  Mathematical  Tables  and 
Other  Aids  to  Computation. 
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4 . Computing  Services 


SEPARATION  OF  EXPONENTIALS 
Project  11.1/31-48-2 


Origin:  MBS 

Sponsor:  Air  Materiel  Command,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  INACTIVE 


TEST  OF  AN  ADAPTATION  OF  THE  MONTE  CARLO  METHOD 
Project  11.1/31-49-10 

Origin:  The  Rand  Corporation  Authorized  11/16/48 

Sponsor:  Air  Materiel  Command,  USAF  Completed  8/31/49 

Managers:  L.  L.  Bailin,  Roselyn  A.  Siegel 

Objective:  By  use  of  an  adaptation  of  the  so-called  "Monte  Carlo" 

method,  to  obtain  a statistical  solution  of  the  neutron  transport  equa- 
tion . 

(1)  'f  + t = | / f (x,  V’)dK', 

and  compare  the  solution  with  one  obtained  by  analytic  methods. 

Background:  One  considerable  disadvantage  in  using  the  Monte  Carlo 

method  to  solve  an  equation  such  as  (l),  where  the  independent  variable 
is  greatly  attenuated  for  certain  physically  interesting  situations,  is 
that  often  an  impossibly  large  number  of  "particles"  must  be  calculated 
in  order  to  obtain  a reasonably  accurate  statistical  approximation  to  the 
answer.  This  project  is  to  test  a conjecture  that  this  difficulty  can 
be  removed  by  substituting  ^ and  solving  the  resulting  equation 

for  / . This  corresponds  to  a physical  situation  where  particles  may 
multiply  as  well  as  be  absorbed. 

Status:  COMPLETED. 

Publication:  Mr.  Herman  Kahn  of  the  Rand  Corporation  expects  to 

publish  his  results  in  Nucleonics . 
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RADIATION  CHARACTERISTICS  OF  A TURNSTILE  ANTENNA 
SHIELDED  BY  A SECTION  OF  A METALLIC  TUBE 
CLOSED  AT  ONE  END 
Project  11.1/31-49-12 

Origin:  North  American  Aviation  Corporation  Authorized  3/^/^9 

Sponsor:  Air  Materiel  Command,  USAF  Completed  9/30/49 

Manager:  L.  L.  Bailin 

Objective:  To  determine  by  theoretical  methods  the  radiation 

characteristics  of  the  antenna,  i.e.,  the  far  field  radiation  pattern 
and  polarization,  the  behavior  of  the  propagating  modes  in  the  cylin- 
drical shielding  can,  the  impedence  seen  by  the  antenna,  the  field  and 
power  gain  functions,  subject  to  the  following  assumptions: 

(i)  The  antenna  consists  of  perpendicularly  crossed  wires, 
of  equal  length,  excited  by  a sinusoidal  current  distribution  which 
differs  90°  in  time  phase  between  members  of  the  cross. 

(ii)  The  shield  consists  of  a cylindrical  can  open  at  one  end. 
The  plane  of  the  antenna  is  perpendicular  to  the  axis  of  the  cylinder; 
the  center  of  the  crossed  wires  is  on  the  axis. 

(iii)  The  shield  is  assumed  to  be  a perfect  conductor;  the 
effect  of  the  antenna  supports  is  neglected. 

Background:  This  problem  was  proposed  by  the  Aerophysics  Labora- 

tory of  the  North  American  Aviation  Corporation  as  a step  in  obtaining  a 
satisfactory  basis  for  efficient  design  of  antennas  of  this  type. 

Comments:  Professors  A.  Banos,  Jr.,  and  D.  S.  Saxon  of  the  Univer- 

sity of  California  at  Los  Angeles  are  acting  as  consultants  on  this 
problem . 

Status:  COMPLETED.  The  results  were  transmitted  to  the  North 

American  Aviation  Corporation. 

Publication:  Paper  by  Mr . L.  L.  Bailin  covering  the  methods  of 

radiation  characteristics  will  be  submitted  for  publication  in  the  NBS 
Journal  of  Research. 


OPTIMUM  STAR  CONFIGURATIONS 
Project  11.1/31-49-1^ 

Origin:  Hughes  Aircraft  Company,  Inc.  Authorized  5/25/^9 

Sponsor:  Air  Materiel  Command,  USAF  Completed  8/8/49 

Manager:  E.  Yowell 

Objective:  To  find  the  parameters  of  certain  types  of  star- 

shaped figures  which  maximize  the  area  between  the  star  and  an  enclos- 
ing circle. 

Background:  This  problem  arises  in  the  study  of  methods  of  pack- 

ing solid  rocket  fuel  in  cylindrical  rockets. 

Status:  COMPLETED.  The  results  were  forwarded  to  the  Hughes 

Aircraft  Corporation.  For  certain  sets  of  parameters  no  solution  exists 
and  mathematical  criteria  were  derived  for  the  determination  of  this  point. 
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COMPUTATIONS  RELATING  TO  THEODOLITE  SYSTEMS 
Project  11.1/31-5P-1 

Origin:  Boeing  Aircraft  Corporation 

Sponsor:  Air  Materiel  Command,  USAF 

Manager:  E.  Yowell 

Objective:  To  compute  certain  quantities  by  formulae  furnished  by 

the  Boeing  .Aircraft  Corporation  for  parameters  corresponding  to  various 
positions  of  a sighted  object. 

Background:  The  problem  arose  in  research  work  of  the  Boeing  Air- 

craft Corporation  under  contract  to  the  Air  Materiel  Command. 

Status:  COMPLETED. 


Authorized  8/1/49 
Completed  8/8/49 


TABLE  RELATING  TO  A LINEAR  ANTENNA  ARRAY 
Project  11.1/31-50-2 

Origin:  Hughes  Aircraft  Corporation 

Sponsor:  Air  Materiel  Command,  USAF 

Manager:  E.  Yowell 

Objective:  To  compute 

f(n,/)  = 10  log10 

for  x = n 7T  /3/36O;  Z = exp  (ijA); 

n = 0(1)135  ; f = 0( 1°) 180°  . 

Background:  The  problem  arose  in  connection  with  a study  of  rela- 

tive currents,  performed  by  Hughes  Aircraft  Corporation  as  part  of  a 
contract  with  the  Air  Materiel  Command. 

Status:  COMPLETED. 


STUDIES  OF  A DIFFERENTIAL  EQUATION 
Project  II.I/3I-5O-3 

Origin:  Hughes  Aircraft  Corporation  Authorized  8/1/49 

Sponsor:  Air  Materiel  Command,  USAF  Completed  9/30/49 

Managers:  Gertrude  Blanch,  E.  Yowell 

Objective:  To  find  the  sinusoidal  solutions  of  the  equation 

y"  - y'  + f2-  = V.  U)  > 

where 

7^  ( x ) =j^sin(sx  + A Z.)  j /xc , A = 0,  1;  c = 0,1;  s = l(l)8;  and  a = 4 . 

To  tabulate  the  function  up  to  a point  where  the  asymptotic  expansion  gives 
approximately  two  significant  figures  . 


1 + e~x 
Z + e~x 


Authorized  8/1/49 
Completed  9/10/49 
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Background:  These  solutions  are  desired  in  connection  with  research 

being  performed  by  the  Hughes  Aircraft  Corporation  under  a contract  with 
the  Air  Material  Command. 

Status:  COMPLETED. 


COMPUTING  SERVICES  FOR  RESEARCH  STAFF 
OF  THE  INSTITUTE  FOR  NUMERICAL  ANALYSIS 
Project  11.1/52-49-1 


Origin:  NBS  Authorized  9/2/48 

Sponsor:  Office  of  Naval  Research,  USN  Revised  11/15/49 

Manager:  Gertrude  Blanch 

Objective:  To  assist  the  research  staff  of  the  Institute  for 

Numerical  Analysis  by  performing  small  computational  tasks. 

Background:  A function  of  the  Computation  Unit  of  the  Institute 

for  Numerical  Analysis  is  to  provide  facilities  for  the  numerical  exper- 
imentation needed  in  the  research  work  of  the  Institute.  Furthermore, 
the  publications  of  the  research  staff  frequently  contain  small  mathe- 
matical tables  and  numerical  demonstrations  of  method,  and  it  is  the  duty 
of  the  Computation  Unit  to  furnish  this  material. 

Status:  CONTINUED.  Among  the  more  extensive  computations  during 

this  quarter  were: 

1.  Application  of  the  method  of  C.  Lanczos  for  the  determina- 
tion of  eigenvalues  and  eigenvectors.  (See  11.1/1-50-5,  section  I). 

2.  Studies  in  the  problem  of  finding  maxima  and  minima  on  sur- 
faces, for  M.  R.  Hestenes . (See  ll.l/l-50-l  and  11.1/1-50-2). 

5.  Studies  in  generating  random  digits,  for  G.  Forsythe.  (See 
11.1/1-49-PM2) 

4.  The  determination  of  "Gibbs’  constant"  for  Holder  means,  for 
0.  Szasz . (See  11.1/1-50-4). 

5.  Approximation  studies  for  A.  M.  Ostrowski . (See  11 . 1/1-49-AE5) . 

6.  Application  of  the  Monte  Carlo  method  to  the  numerical  sol- 

ution of  the  Dirichlet  problem,  for  harmonic  functions  in  a rectangle; 
for  W.  Feller:  From  any  point  in  the  boundary  of  a region  a random  walk 

is  made  along  a square  grid.  The  walk  is  stopped  when  a boundary  is 
reached,  and  the  value  of  the  function  at  the  boundary  point  is  computed. 
The  average  value  of  the  boundary  function  for  a large  number  of  walks 
gives  the  value  of  the  function  at  the  starting  point. 

7.  Obtaining  the  limit  An,  as  n approaches  infinity,  for  a stoch- 
astic matrix  A by  calculating  the  projective  operator  of  the  space, 
complementary  to  that  spanned  by  the  vectors  of  the  matrix  (A  - I) . The 
problem  was  performed  by  H.  Weinberger,  under  the  direction  of  W.  Feller. 

8.  Computations  relating  to  the  zeros  of  the  L-series,  for 
J.  Barkley  Rosser. 

9.  Problems  in  number  theory,  performed  by  Robert  Douthitt  at 
the  request  of  D.  H.  Lehmer  on  IBM  equipment.  Mr.  Douthitt  punched  and 
accumulated  the  quotients  of  the  first  10,000  integers.  He  likewise 
punched  and  accumulated  the  values  of  5n^  for  the  same  range.  A predic- 

2 

tion  of  D.  H.  Lehmer  was  not  borne  out77^  by  the  experiment. 

10.  Study  of  round-off  error,  for  G.  Forsythe  and  E.  Yowell.  Step- 
wise integrations,  of  considerable  length,  of  the  sine  and  cosine  func- 
tions were  performed  using  a certain  method  of  random  round-off . 
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11.  Computation  of  lowest  eigenvalue  of  Schrodinger ' s equa- 
tion by  a Monte  Carlo  method,  for  M.  Kac,  M.  Donsker  and  E.  Yowell. 

See  project  11.1/1-50-3,  section  (IV). 


REDUCTION  OF  THEODOLITE  DATA 
Project  11.1/32-49-2 

Origin:  Naval  Air  Missile  Test  Center  (Point  Mugu) 

Sponsor:  Bureau  of  Aeronautics,  USN 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  INACTIVE. 


THE  DETERMINATION  OF  THE  PERIODS  AND  AMPLITUDES 
OF  THE  LIGHT  VARIATIONS  OF  THE  STARS 
& SCUTI  AND  12  LACERTAE 
Project  11.1/32-49-4 


Origin:  MBS 

Sponsor:  Office  of  Air  Research,  AMC , USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  INACTIVE 


DETERMINATION  OF  ORBITS  OF  COMETS, 
MINOR  PLANETS,  AND  SATELLITES 
Project  11.1/32-49-6 


Origins  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  About  35$  completed.  Work  was  continued  by 

Dr.  S.  Herrick  and  his  staff  at  the  University  of  California  at  Los 
Angeles,  using  Institute  for  Numerical  Analysis  I.B.M.  equipment,  when 
machine  time  was  available. 


Status  of  Projects 
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TEST  OF  A VARIATIONAL  PROCEDURE  FOR  DETERMINING  THE  LOWEST  ENERGY 
VALUE  OF  A SIMPLE  QUANTUM  MODEL 
Project  11 .1/32-50-1 


Origin:  NBS  Authorized  8/1/49 

Sponsor:  Office  of  Naval  Research,  USN 

Manager:  H.  Gruen 

Objective:  To  test  the  feasibility  of  a variational  method  for 

computing  quantum  states  by  applying  the  method  to  the  following  problem 
in  nuclear  physics. 


The  operator 


H 


h2  d2 
2m  dx 


2 + V0 


hid£ 

2m  dx2 


M - a , 

|x|  ^ a , 


is  the  Hamiltonian  for  a particle  of  a mass  m in  a potential  well  of  a 
depth  Vq  and  a width  2a.  h is  Planck's  constant  divided  by  2tt  . 

A trial  wave  function  is  chosen  as 


$(x) 


9 


|x|  ^ a , 
|x|  - a . 


We  define 

4n<A> 

Fn<A>  = A + TA)  ’ 

• n -1 


where  n = 1,3,5,...,  and 

/3n(A)  = / § (x)  (H  - A) n i (x)dx//  $(x)j  2 dx. 

The  minimum  values  of  F(A)  for  A ranging  from  -500  to  +500  for 
odd  values  of  n 4 9 will  be  computed  and  the  results  compared  with  the 
theoretical  conjecture  (see  Background)  that,  an  optimum  n0  exists  where 
mj^n  FnQ(A)  is  the  best  approximation  to  E0  obtainable  with  the  given  trial 
wave  function. 


Background:  The  variational  method  together  with  the  conjectured 

results  was  suggested  by  Dr.  Edward  Gerjuoy  of  the  University  of  Southern 
California . 


Status:  NEW 
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CHARACTERISTIC  VALUES  OF  MATHIEU'S  DIFFERENTIAL  EQUATION 

Project  11.1/4-49-6 


Origin:  NBS 

Sponsor:  Air  Materiel  Command,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Characteristic  values  were  being  computed  for 

s beyond  the  range  of  the  published  tables,  as  outlined  in  previous 
status  reports.  At  the  end  of  the  quarter,  be^(s)  was  available  at  an 
interpolable  interval  in  1/ N/s,  from  0 to  .1  (i.e.  for  s from  100  to  oo ) . 


SUMMER  FELLOWSHIP  EDUCATIONAL  PROGRAM 
Project  11.1/4-50-1 

Origin:  NBS  Authorized  7/1/49 

Sponsor:  Office  of  Naval  Research,  USN  Completed  8/26/49 

Managers:  E.  C.  Yowell,  H.  D.  Huskey,  W.  Feller,  A.  M.  Ostrowski 

Objective:  To  acquaint  a selected  group  of  graduate  students  with 

some  of  the  trends  and  techniques  of  modern  numerical  analysis. 

Background:  The  Summer  Fellowship  program  was  suggested  by  Pro- 

fessors Morgan  Ward  and  H.  F.  Bohnenblust  of  the  California  Institute  of 
Technology . 

Status:  COMPLETED.  The  following  courses  were  given: 

(l)  "Techniques  of  present-day  computing".  Application  of 
electronic  punched  card  equipment  to  solution  of  mathematical  problems. 

Dr.  E.  C.  Yowell  and  associates.  (2)  "Theory  of  automatic  digital  computing 
machinery".  Logical  aspects  of  the  design  of  present  and  proposed  systems; 
programming  and  problem  preparation.  Dr.  H.  D.  Huskey  and  associates. 

(5)  "Seminar  in  probability  methods."  Introduction  to  the  partial  dif- 
ferential equations  which  arise  in  the  theory  of  stochastic  processes,  and 
study  of  sampling  methods  for  solving  these  equations.  Dr.  Feller  and 
members  of  the  staff.  (4)  "Seminar  in  the  Inversion  of  matrices  of  high 
order."  Chiefly  a study  of  the  paper,  "Numerical  inverting  of  matrices 
of  high  order,"  by  J.  von  Neumann  and  H.  H.  Goldstine,  Bui.  Amer . Math. 

Soc.,  52,  1021  to  1099,  (1947).  Dr-  A.  M.  Ostrowski. 

The  summer  fellows  were:  Robert  J.  Diamond,  California  Insti- 
tute of  Technology;  Robert  C.  Douthitt,  University  of  California,  Berkeley; 
Ernest  Elyash,  Cornell  University;  Harold  Gruen,  University  of  California, 

Los  Angeles;  James  G.  C.  Templeton,  Princeton  University;  and  Hans  Weinberger, 
Carnegie  Institute  of  Technology. 


II.  Computation  Laboratory 
(Section  11.2) 


1 . Mathematical  Tables 


TABLES  OF  E1(z),  (z  = x + iy) 
Project  11.2/2-43-3 


Origin:  Canadian  National  Research  Council 

Sponsor:  Office  of  Air  Research,  AMC , USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  completed;  preparation  of  final 

manuscript  for  publication  was  continued. 


TABLE  OF  THE  GAMMA  FUNCTION  FOR  COMPLEX  ARGUMENTS 
Project  11.2/2-46-1 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  INACTIVE 


MATHIEU  FUNCTIONS,  I 
Project  11.2/2-46-2 


Origin:  Applied  Mathematics  Panel,  NDRC  Authorized  7/1/47 

Manager:  Gertrude  Blanch  Completed  9/15/^ 9 

Objective:  To  compute  an  eight-place  table  of  the  first  13  odd  and 

16  even  characteristic  values  b of  Mathieu ' s differential  equation 

y"  + (b  - s cos^t)y  = 0 

for  s ranging  from  0 to  100  at  various  intervals,  and  the  Fourier  coeffi- 
cients of  the  solutions  corresponding  to  these  characteristic  values,  as 
well  as  certain  related  functions. 
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Background:  Mathieu  functions  arise  in  the  solution  of  the  wave 

equation  for  elliptical  domains . Numerous  physical  applications  involving 
Mathieu  functions  are  described  in  "Theory  and  Applications  of  Mathieu 
functions,"  by  N.  W.  McLachlan  (Oxford  Press,  1947).  The  existing  tables 
by  Stratton-Morse-Chu-Hutner,  Ince  and  Goldstein  are  inadequate.  The 
project  was  proposed  by  Dr.  Philip  Morse  and  Dr.  James  Wakelin. 

Comments:  Related  to  project  11.1/2-45-1  ('formerly  11.2/2-45-1,  prior 

to  April  1949) • 

Status:  COMPLETED. 

Publication:  Manuscript  transmitted  to  Columbia  University  Press 

for  publication. 


TABLES  OF  COULOMB  WAVE  FUNCTIONS 
Project  11.2/2-47-2 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC , USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  completed';  checking  was  continued. 


TABLE  OF  ANTILOGARITKMS 
Project  11.2/2-47-5 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  completed;  checking  was  continued. 


TABLES  FOR  THE  OCCASIONAL  COMPUTER 
Project  11.2/2-47-4 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 


Status:  INACTIVE 
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TABLE  OF  LAGRANGIAN  COEFFICIENTS 
FOR  SEXAGESIMAL  INTERPOLATION 
Project  11.2/2-48-2 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  INACTIVE 


ZEROS  AND  WEIGHT  FACTORS  OF  THE  FIRST  SIXTEEN 
HERMITE  POLYNOMIALS 
Project  11.2/2-49-1 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  INACTIVE 


RADIX  TABLE  FOR  CALCULATING  LOGARITHMS  TO  MANY  PLACES 
Project  11.2/2-49-2 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue 

Status:  INACTIVE 


TABLE  OF  POWERS  OF  COMPLEX  NUMBERS 
Project  11.2/2-50-1 


Origin:  NBS 

Sponsor:  Office  of  Air  Research,  AMC,  USAF  Authorized  9/l/^9 

Manager:  Irene  Stegun 

Objective:  To  publish  a table  of  the  exact  values  of  the  first 

25  powers  of  z = x + iy  for  x,  y = 0(1)10. 

Background:  The  manuscript  of  such  a table,  which  is  basic  for  a 

great  many  computations,  had  been  prepared  several  years  ago  by  the 
Mathematical  Tables  Project,  under  Mr.  Herbert  Salzer.  At  the  present 
time  it  is  intended  to  check  and  publish  this  manuscript. 

Status:  NEW.  Checking  of  galley  proofs  was  in  progress. 

Publication:  The  table  Is  being  printed  by  the  Government  Printing 
Office  and  will  be  issued  as  No.  8 in  the  NBS  Applied  Mathematics  Series. 
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5* 


TABLES  TO  FACILITATE  SEQUENTIAL  TESTS 
Project  11.2/2-50-2 


Origin;  NBS  Authorized  9/l/^9 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Manager:  Irene  Stegun 


Objective : 
equation : 


for  n = 1(1)200; 


To  publish  a table  of  the  solution 

2 


L = log  F 


[n  1. 
[2’  2 ’ 


S2z1 


nS 

2 


S = .1( .1)1(  .2)2,  2.5;  ± L = 2(1)7, 


z(L,n,<5)  of  the 


In  19,  In  99. 


Background:  These  tables  were  originally  computed  by  the  Mathematical 

Tables  Project  for  the  Statistical  Research  Group  at  Columbia  University 
during  the  war  for  use  in  connection  with  sequential  testing.  They  had  not 
been  published  before. 


Status:  NEW.  Checking  of  galley  proofs  was  in  progress. 

Publication:  The  table  is  being  printed  by  the  Government  Printing 
Office  and  will  be  Issued  as  No.  7 in  the  NBS  Applied  Mathematics  Series. 


TABLE  OF  CHEBYSHEV  POLYNOMIALS 
Project  11.2/2-50-5 


Origin:  NBS  Authorized  9/l/^9 

Sponsor:  Office  of  Air  Research,  AMC,  USAF 

Manager:  Irene  Stegun 

Objective:  To  publish  a table  of  the  Chebyshev  polynomials 

Sn(x),  Cn(x)  for  n = 2(1)12  and  x = 0(. 001)2. 

Background:  The  manuscript  of  this  table  was  prepared  some  time  ago 

by  the  Mathematical  Tables  Project  but  publication  was  postponed.  These 
polynomials  are  frequently  used  as  approximations  to  more  complicated 
functions.  They  are  especially  convenient  for  use  in  connection  with 
high-speed  automatic  computing  machines. 

Status:  NEW.  Checking  of  galley  proofs  was  in  progress. 

Publication:  The  table  is  being  printed  by  the  Government  Printing 
Office  and  will  be  issued  a3  No.  9 in  the  NBS  Applied  Mathematics  Series. 
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2 . Computing  Services 


RANGE  LAG  OF  BOMBS 
Project  11.2/51-49-2 


Origin:  Air  Materiel  Command,  USAF  Authorized  6/25/49 

Sponsor:  Office  of  Air  Research,  AMC,  USAF  Completed  8/1/49 

Manager:  F.  L.  Alt 

Objective:  To  approximate  the  range  lag  of  bombs  as  a simple 

analytic  function  of  altitude  and  air-speed  at  time  of  release,  and  of 
ballistic  coefficient. 

Background:  The  range  lag,  i.e.,  the  difference  between  the  actual 
range  (horizontal  travel)  of  a bomb  and  that  range  which  would  occur  in  a 
vacuum.  Is  fitted  by  means  of  least  squares  by  an  expression  of  form 
My  + Ny?/2,  where  y is  the  altitude  of  release  and  where  M and  N are  co- 
efficients which  in  turn  are  to  be  fitted  by  simple  functions  of  air-speed 
and  ballistic  coefficient.  The  reports  are  desired  by  contractors  of 
the  Air  Materiel  Command  for  the  design  of  certain  compensating  equipment. 

Status:  COMPLETED.  Manuscript  transmitted  to  the  Air  Materiel 

Command . 


GUST  ATTACKS  ON  DELTA  WING 
Project  11.2/31-50-1 

Origin:  Aircraft  Laboratory,  AMC,  USAF  Authorized  9/15/49 

Sponsor:  Office  of  Air  Research,  AMC,  USAF  Estimated  Completion  12/49 

Manager:  Irene  Stegun 

Objective:  To  evaluate  by  numerical  computation  the  basic  un- 
steady variations  in  lift  caused  by  a sharp-edged  vertical  wind  gust  at- 
tacking a portion  of  a delta  wing  in  supersonic  flight. 

Status:  NEW.  Preliminary  computations  were  started. 


HEAT  CONDUCTION  EQUATION 
Project  11.2/33-46-1 


Origin:  Bureau  of  Ordnance,  Department  of  the  Navy 

Full  project  description  appears  in  Apr-Jun  1949  issue. 


Status:  INACTIVE 
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FOURIER  TRANSFORM  ADJUSTMENT  COMPUTATIONS 
Project  11.2/33-49-2 


Origin:  Naval  Research  Laboratory 

Sponsor:  Naval  Research  Laboratory 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  for  various  values  of  the  para- 

meter ap,  mp,  were  performed  when  requested. 


TABLES  OF  THERMODYNAMIC  PROPERTIES  OF  GASES 
Project  11.2/33-49-5 


Origin:  Section  3*2,  NBS 

Sponsor:  National  Advisory  Committee  on  Aeronautics 

Full  project  description  appears  In  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Tables  of  zero  pressure  properties  were  com- 

pleted for  atomic  oxygen,  atomic  hydrogen,  and  atomic  nitrogen.  Tables 
of  free  energy  properties  were  completed  for  No,  Op,  normal  Hp,  atomic 
oxygen,  atomic  nitrogen  and  atomic  hydrogen.  A table  of  vapor  pressure 
of  oxygen  was  also  completed.  Prior  to  this  quarter,  tables  of  zero 
pressure  properties  had  been  completed  for  Np,  Op,  A,  air,  and  normal 
Hp;  and  computations  for  real  gas  properties  had  been  completed  for  dry 
air,  normal  Hp,  and  moist  air. 


ADMINISTRATIVE  RECORDS  ON  PUNCHED  CARDS 
Project  11.2/33-49-8 


Origin:  Division  M,  NBS  Authorized  12/15/48 

Sponsor:  Division  M,  NBS  Terminated  9/15/49 

Manager:  A.  H.  Rosenthal 

Objective:  To  establish  and  maintain  punched-card  records  of  cer- 

tain administrative  items,  such  as  stock  records,  travel  data,  etc.  To 
prepare  up-to-date  printed  listings  of  the  records  as  required. 

Background:  Such  records  were  formerly  kept  by  manual  methods. 

The  punched  card  method  will  make  it  possible  to  compile  necessary  per- 
iodic reports  in  much  less  time  and  will  make  record  data  more  easily 
and  quickly  available  for  use  in  statistical  studies. 

Status:  TERMINATED.  The  project  was  transferred  to  the  NBS  Budget 

and  Management  Division  (M) . 
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FERMI  FUNCTION  II 
Project  11.2/33-49-10 


Origin:  Section  4.4,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  completed  and  checked.  Prepara- 

tion of  auxiliary  tables  and  of  an  introduction  to  accompany  publication 
of  the  table  was  started. 

Publication:  To  be  submitted  for  publication  in  the  NBS  Applied 

Mathematics  Series . 


EQUILIBRIUM  COMPOSITION  OF  COMBUSTION  GASES 

Project  11.2/33-49-11 


Origin:  Lewis  Flight  Propulsion  Laboratory,  NACA 

Financial  Sponsor:  Lewis  Flight  Propulsion  Laboratory,  NACA 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  were  started. 


SHOCK  WAVE  PARAMETERS 
Project  11.2/33-49-13 


Origin:  Bureau  of  Ordnance,  Department  of  the  Navy 

Sponsor:  Bureau  of  Ordnance,  Department  of  the  Navy 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  were  completed,  and  preparation 

of  the  manuscript  for  publication  of  the  tables  of  In(x,cv)  was  under  way. 
An  investigation  of  the  interpolability  of  these  tables  was  initiated. 


BASIC  IONOSPHERIC  DATA 
Project  11.2/33-49-14 


Origin:  Section  14.1,  NBS 

Sponsor:  Section  14.1,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Recording  and  processing  of  the  observation 

data  was  continued. 
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Publications:  Current  data  are  published  in  the  monthly  CRPL-F 

reports  of  the  NBS  Central  Radio  Propagation  Laboratory;  copy  prepared 
from  punched  cards  will  replace  hand-made  copy  as  soon  as  procedures 
and  schedules  are  worked  out.  The  results  of  statistical  computations 
will  be  included  in  scientific  papers  where  appropriate . 


SEQUENTIAL  SAMPLING  PLANS 
Project  11.2/33-49-15 


Origin:  Hq . Quartermaster  Corps  Inspection  Service, 

Department  of  the  Army 

Sponsor:  National  Applied  Mathematics  Laboratories 

Manager:  Irene  Stegun 

Objective:  To  design  truncated  sequential  sampling  plans  which  have 

OC  curves  matching  those  of  a given  set  of  single  sampling  plans  at  two 
points . 

Background:  This  set  of  sequential  sampling  plans  is  to  be  used  in 

the  revision  of  JAN  Standard  105. 

Comments:  Supervision  is  carried  out  in  close  contact  with 

Mr.  J.  Cameron  of  the  NBS  Statistical  Engineering  Laboratory. 

Status:  COMPLETED.  Manuscript  transmitted  to  the  originator. 


Authorized  4/14/49 
Completed  9/30/49 


STRAIN  ROSETTE  MEASUREMENTS 
Project  11.2/33-49-16 


Origin:  David  Taylor  Model  Basin,  USN  Authorized  4/14/49 

Sponsor:  David  Taylor  Model  Basin,  USN  Completed  9/15/49 

Manager:  Irene  Stegun 

Objective:  To  evaluate  a series  of  strain  rosette  measurements 

taken  at  the  David  Taylor  Model  Basin. 

Background:  "Strain  rosette  measurements"  consist  in  determining 

strains  in  three  coplanar  directions,  two  being  perpendicular  to  each 
other  and  the  third  at  a 45°  angle.  From  these  measurements  it  is  de- 
sired to  determine  the  magnitude  and  directions  of  the  principal  strains. 

Status:  COMPLETED.  Manuscript  transmitted  to  the  David  Taylor 

Model  Basin. 


Status  of  Projects 


37 


RADIO -TELEGPAPH  INTERFERENCE 
Project  11.2/33-49-17 

Origin:  Section  14.4.  NBS 

Sponsor:  Section  14.4,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Recording  of  data  was  continued. 


PERCENTAGE  POINTS  OF  THE  ARITHMETIC  MEAN  IN  RANDOM  SAMPLES 
FROM  THE  SECH  AND  SECH2  DISTRIBUTIONS 
Project  11.2/33-49-18 


Origin:  Section  11. 3,  NBS 

Sponsor:  Section  11.3,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  about  75$  completed. 


COMPUTATION  OF  NORMAL  MODES  OF  BOX  BEAM 
Project  11.2/33-49-20 


Origin:  Section  6.4,  NBS 

Sponsor:  Section  6.4,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  CONTINUED.  75$  completed. 


TABLES  OF  CABLE  FUNCTIONS 
Project  11.2/33-^9-21 


Origin:  David  Taylor  Model  Basin,  USN 

Sponsor:  David  Taylor  Model  Basin,  USN 

Manager:  S.  Prusch 


Authorized  6/25/49 
Completed  8/49 


Objective : 


To  evaluate  the  functions 

loge  r = J ^ d0  > 

o ^ 


r c os/h/ 


‘ =/ 
D 
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where  P 


f — CO30  + w sin  0,Q,=  - /sin/ I sin/  + 
/ COs/|  9 ' 


w CO30  , 


and  w = sec  0C  - cos  0C , for  f = .0l(  .01)  .05,  jZfc  = 0(5°)90°, 

/ = [/c(l°)/c  +180°]. 

Background:  These  functions  facilitate  the  determination  of  the 

shape  of  and  tensions  in  a flexible  cable  that  is  held  in  a uniform 
stream.  The  effects  of  the  weight  of  the  cable  and  frictional  drag  are 
included.  The  types  of  configurations  that  may  be  treated  by  means  of 
these  functions  are  most  general. 

Status:  COMPLETED.  Manuscript  transmitted  to  the  David  Taylor 

Model  Basin. 


WAVE  RESISTANCE  OF  SHIPS 
Project  11.2/55-50-1 


Origin:  David  Taylor  Model  Basin,  USN  Authorized  9/l/^9 

Sponsor:  David  Taylor  Model  Basin,  USN 

Manager:  J.  H.  Levin 

/ YL 

Objective:  To  calculate  functions  Mn(y)  = J £ x^  sink$  d!, 

Mn'(k)  = Jo  f cos  V | dj  for  n = 0(  1)5(2)  11,  and  for  approximately  150 
values  of  V in  the  range  l^Y^lOO;  to  calculate  also  specified  products 
of  the  type  MpM^,  Mj_  'Mk'.  The  above  integrals  are  evaluated  by  the  use 
of  the  recursion  formulae: 


M (K)  = _ 003  y + n alnV  _ n(n-l)  Mn_ 


v 

M <(y)  = 3ln  y + n 003  V _ P.tP-1!  M'  o 

n \ ) y y 2.  yz  n -2 

Background:  For  simple  forms  the  wave  resistance  of  a ship  is  given 

to  a first  approximation  by  Michell's  integral: 


r = 8^  bV  /{ i+2(y)  +^20')}v/tm7F_i 

o 

' r1 

I+(V)  = J L n an$  n-X  cos  Y J d i J (l  -$  m)  e'0*  dt 

o 0 

' I 

I+(  k)  = _/  21  n an  $ n-l  sin  V ! d ^ J (l  - t ™)  e‘0i 


(1) 

(2) 

(5) 


y/y*v 


where  n,  m are  integers,  the  quantity  preceding  the  integral  sign  in  (l)  is 


a dimensional  factor,  and  9 = — - — ( ? on 3 tan t ) . This  approximation 


Vo 


can 


be  used  as  a basis  for  more  complicated  cases  . An  equivalent  expression 
(Havelock's  Integral)  is: 
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A 2 r,3( 

R = S.  J jl+ 2 ( Y0  sec  9)  + Y0  sec  9)|  sec^Q  d9 

where  Y0  sec  9 Is  substituted  for  V . The  expressions  Mn(k)  and  Mn'(y)  are 
necessary  for  the  evaluation  of  either  Michell's  integral  or  Havelock's 
integral . 

The  indices  n of  the  functions  Mn(y),  Mn'(K)  depend  upon  expon- 
ents of  terms  in  polynomials  used  for  approximating  ship  forms.  With  the 
chosen  indices  (exponents)  the  whole  range  of  normal  ship  forms  can  be 
covered. 

Besides,  the  same  functions  Mn(y),  Mn'(y)  admit  of  solving  resistance 
problems  of  totally  submerged  bodies  near  the  surface,  gliding  vessels,  and 
problems  of  behavior  of  a ship  in  a.  seaway . 

Thus  these  functions  are  of  basic  importance  in  Theoretical  Naval 
Architecture . 

Status:  NEW.  Computations  were  started. 


RATING  OF  WATER  CURRENT  METERS 
Project  11.2/33-50-2 


Origin:  Division  6.5,  NBS  Authorized  9/15/49 

Sponsor:  Division  6.5,  NBS 

Manager:  J.  H.  Levin 

Objective:  To  represent  the  relationship  between  velocity  of  current 

and  meter  reading  by  a fitted  curve,  which  is  to  be  used  in  calibrating 
the  meter. 

Background:  In  rating  current  meters,  a series  of  towing  tests  is 

made  in  still  water.  The  meter  is  towed  through  the  water  at  a known 
velocity  and  the  revolutions  per  second  of  the  meter  recorded.  The  curve 
representing  the  velocity  is  called  the  rating  curve.  For  purposes  of 
calibration  the  rating  curve  may  be  approximated  by  two  intersecting  straight 
lines.  The  intersection  occurs  at  about  the  same  meter  reading  for  each 
type  of  meter.  These  lines  are  used  in  the  calibration  of  the  meters..  The 
equations  of  the  lines  are  obtained  from  the  .observed  data  on  punched  card 
machines.  Also  obtained  are  the  deviations  of  the  observed  current  veloci- 
ties from  the  current  velocities  calculated  from  the  lines  . 

Status:  NEW.  Computations  were  performed  as  requested. 
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A PROBLEM  IN  MOLECULAR  STRUCTURE 
Project  II.2/33-50-3 


Origin:  Naval  Research  Laboratory,  USN 
Sponsor:  Naval  Research  Laboratory,  USN 
Manager:  Irene  Stegun 

Objective:  To  tabulate  values  of  the  function 


Authorized  9/1/49 


f = e 


for  q > 0, 


f 

vhere 


e x/lqTi  [I_i.  (f“)  + I 

= e,ah+1<R"L>  [l_l  (§V  U l3!)].  for 


q = 


L h ( l-4a  /3)  - R (h+/3  ) 


(4  ah+l)2  [h(4ah+l)  - 4 a (h+/3)"| 


a and  L are  assigned  constants,  and  R ranges  from  3.9  to  4.6,  /3  from  1 
to  7,  and  h from  20  to  100. 

Background:  This  function  occurs  in  investigations  of  molecular 

structure.  It  represents  the  probability  distribution  of  the  distance  R 
of  two  atoms  in  the  internal  vibratory  motion  of  the  molecule,  where  L 
is  the  equilibrium  distance  of  the  atoms  at  rest  and  h,  /3  are  damping 
factors  relating  to  the  vibratory  motion. 

Comments:  This  problem  was  formulated  by  Dr.  J.  Karle  and 

Mr.  H.  Hauptman  of  the  Naval  Research  Laboratory. 

Status:  NEW 


STANDARD  L0RAN  TABLES 
Project  11.2/34-50-1:  Gulf  Coast  Chain 


Origin:  U.  S.  Navy  Hydrographic  Office  Authorized  9/1/49 

Sponsor:  U.  S.  Navy  Hydrographic  Office 

Manager:  J.  H.  Levin 

Objective:  Preparation  of  tables  giving  coordinates  of  hyperbolic 

lines  of  positions. 

Background:  Standard  Loran  Navigation  tables  are  necessary  for 

preparation  of  charts  used  by  navigators  in  determining  their  positions 
with  the  aid  of  certain  electronic  equipment. 

Status:  NEW.  Preliminary  computations  started. 
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EFFECT  OF  NUCLEAR  RADIATIONS  ON  HUMAN  BEINGS 
Project  11.2/55-49-1 


Origin:  Operational  Research  Office,  U.S.Army  (John  Hopkins  University) 

Sponsor:  Operational  Research  Office,  U.S.Army  (John  Hopkins  University) 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  of  attenuation  of  radiation  for 

additional  energies  were  started. 


TABLES  FOR  COMBAT  CREW  PLANNING 
Project  11.2/56-49-2 


Origin:  Air  Comptroller's  Office,  USAF 

Sponsor:  Air  Comptroller's  Office,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computations  had  been  completed  and  checked  on 

punched  cards.  Preparation  of  final  manuscript  continued. 

Publication:  Tables  to  be  included  by  the  Air  Comptroller's  Office 

in  the  next  edition  of  the  Air  Force  Technical  Manual  on  Combat  Crew 
Planning . 


A PROBLEM  IN  LINEAR  PROGRAMING 

Project  11.2/56-49-5 


Origin:  Air  Comptroller's  Office,  USAF 

Sponsor:  Air  Comptroller's  Office,  USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 


Status.  CONTINUED.  Several  satisfactory  procedures  for  com- 
pletely solving  this  type  of  problem  on  IBM  machines  have  been  devised 
and  rested.  Present  efforts  are  being  directed  to  the  attainment  of 
greater  speed  of  solution  and  greater  generality  in  the  type  of  problem 
handled . 


MAGNETIC  MOMENT  OF  ELECTRONS 
Project  11.2/57-49-1 

Origin:  Atomic  Energy  Commission  ' Authorized  6/50/49 

Sponsor:  Atomic  Energy  Commission  Completed  9/15/49 

Manager:  M.  Oliphant 

Objective:  To  evaluate  analytically  a certain 

tegral  of  a rational  function  of  four  variables. 


four-fold  in- 
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Background:  The  integral  represents  effects  of  fourth-order  elec- 

tromagnetic interactions  on  the  magnetic  moment  of  electrons. 

Comments:  The  problem  was  proposed  by  Drs . R.  Karplus  and  M.  Kroll 

of  the  Institute  for  Advanced  Study,  Princeton. 

Status:  Terminated  - at  the  request  of  the  originators,  who  are 

reconsidering  the  scope  of  the  problem. 


3 . Training 


TRAINING  OF  PERSONNEL  IN  THE  OPERATION  OF  NEW  IBM  MACHINES 

Project  11.2/4-49-3 


Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  CONTINUED. 


III.  Statistical  Engineering  Laboratory 
(Section  11.3) 


1.  Research  in  Mathematical  Statistics 


THE  MEAN  DEVIATION,  STANDARD  DEVIATION,  AND  RANGE  AS 
ESTIMATORS  OF  SCALE  PARAMETERS  (MEASURES  OF 
DISPERSION)  OF  PROBABILITY  DISTRIBUTIONS 
Project  11.3/1-47-2 


Origin:  NBS 

Full  project  description  appears  in  Oct-Dec  1948  issue. 

Status:  INACTIVE. 

Publications:  (l)  Probability  center  lines  for  standard  devia- 

tion and  range  charts.  Churchill  Eisenhart.  Industrial  Quality  Con- 
trol VI,  No.  1,  24-26  (July  1949).  (2)  The  relative  frequencies  with 

which  certain  estimators  of  a normal  population  tend  to  underestimate 
its  value.  Churchill  Eisenhart  and  Celia  S.  Martin.  To  be  submitted 
to  the  Journal  of  the  American  Statistical  Association.  An  abstract 
appears  in  Annals  of  Mathematical  Statistics  XIX,  No.  1,  600  (Dec. 1948). 


STATISTICAL  TESTS  OF  SIGNIFICANCE  FOR 
WHEN  THE  NUMBER  OF  OBSERVATIONS 
Project  11.3/1-47-3 


2x2  TABLES 
IS  SMALL 


Origin:  NBS 

Full  project  description  appears  in  Oct-Dec  1948  issue. 


Status:  INACTIVE. 


STATISTICAL  PROPERTIES  OF  SAMPLES  OF  THREE  OBSERVATIONS 

Project  11.3/1-49-1 


Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  About  90$  completed.  Preparation  of  Pub- 

lication (l),  giving  a complete  discussion  of  the  theoretical  and 
numerical  results  developed  under  this  pro ject, approached  completion. 

Publications:  (l)  Properties  of  statistics  involving  the  closest 

pair  in  a sample  of  three  observations.  J.  Lieblein.  To  be  submitted 
to  the  NBS  Journal  of  Research.  (2)  The  fallacy  of  the  best  two  out 
of  three.  For  this  summary  of  the  results  of  the  study  and  their  prac- 
tical implications  see  Tech.  News  Bull.  35,  No.  8,(1  July  1949). 
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STATISTICAL  PROCEDURES  FOR  INTERPOLATED  MEDIANS 
Project  11.3/1-49-2 

Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  INACTIVE 


ELEMENTARY  THEORY  OF  STOCHASTIC  PROCESSES 
Project  11.3/1-49-3 

Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Manager:  Henry  B.  Mann.  (Dr.  Mann, who  is  responsible  for  technical 

aspects  of  this  project,  undertook  the  project  during  a temporary 
employment  by  the  Bureau  from  March  to  June  1949,  and  is  expected  to 
return  for  brief  periods  until  the  work  is  finished.  During  his  absence 
Dr.  Eisenhart  handles  inquiries  and  other  matters  related  to  the  project.) 

Status:  CONTINUED.  A first  draft  of  a monograph  on  elementary 

theory  of  stochastic  processes  had  been  prepared  by  Dr.  Mann,  assisted 
by  Dr.  Miriam  L.  Yevick.  It  consists  of  six  chapters : 1)  Stochastic 

processes  with  a continuous  parameter,  2)  The  Markoff  process  with  con- 
tinuous increment,  3)  The  Fourier  analysis  of  stochastic  processes, 

4)  Stationary  Markoff  processes,  5)  Markoff  processes  with  discontinuous 
Increment,  and  6)  Differential  processes  modified  by  mechanical  devices. 

During  July,  a limited  number  of  copies  of  the  draft  was  cir- 
culated to  certain  selected  individuals  for  their  comments.  It  is  anti- 
cipated that  Dr.  Mann  will  return  in  December  1949  to  revise  the  document 
in  the  light  of  comments  received  and  to  consider  in  detail  the  inclusion 
of  references,  illustrative  examples,  and  Instructional  aids.  Dr.  Mann 
plans  to  use  the  monograph  in  a course  on  stochastic  processes  that  he  will 
give  next  spring  as  a Visiting  Professor,  at  the  University  of  California, 
Berkeley.  It  is  anticipated  that  a final  draft  of  the  monograph  will  then 
be  prepared . 


2.  Manuals  of  Statistical  Methods 


FORMULAS  FOR  OPERATING  CHARACTERISTICS  AND  SAMPLE  SIZES 
FOR  CERTAIN  STATISTICAL  TESTS 
Project  11.3/2-47-2 

Origin:  NBS 

Full  project  description  appears  in  Jan -Mar  1949  issue. 


Status : 


INACTIVE. 
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STANDARD  SAMPLING -INSPECTION  PROCEDURES 
Project  11.3/2-48-1 

Origin:  Office  of  Naval  Research,  and  Research  and  Development  Division 

of  the  Department  of  the  Army 
Full  project  description  appears  in  Apr-Jun  1949  issue. 

Comments:  In  order  to  foster  a wider  appreciation  and  understanding 

of  some  of  the  standard  sampling -inspection  systems  now  in  use  in  the 
Department  of  Defense  and  elsewhere,  a one -term,  in-hours  course  in  the 
basic  statistical  principles,  devices,  and  techniques  involved  was  offered 
at  the  Bureau  (project  11.3/4-49-7)  during  the  period  Mar-Jun  1949.  See 
Projects  and  Publications  Apr-Jun  1949* 

Status:  CONTINUED.  Work  on  this  project  from  April  1949  to  date 

took  the  form  of  active  participation  by  Dr.  Eisenhart  in  the  work  of  the 
Sub-Committee  for  Revision  of  JAN  Standard  105,  which  is  a sub-committee 
of  the  Sampling  Inspection  Committee  of  the  Munitions  Board  Material  In- 
spection Agency.  The  full  membership  of  this  Sub-Committee  consists  of 
Mr.  David  H.  Schwartz  (Quartermaster  Corps  Inspection  Service),  chairman, 
Churchill  Eisenhart  (National  Bureau  of  Standards),  Lee  Goulden  (inspection 
Administration,  Navy  Department),  0.  C.  Griffith  (Air  Materiel  Command), 

Frank  E.  Grubbs  (Ballistics  Research  Laboratories,  Aberdeen  Proving  Ground), 
Mary  J.  Natrella  (Bureau  of  Ships,  Navy  Department),  V.  R.  Pabst,  (Bureau 
of  Ordnance,  Navy  Department),  and  E.  W.  Rote  (Air  Materiel  Command).  Each 
member  serves  solely  as  an  expert  in  his  field,  not  as  an  official  repre- 
sentative of  his  agency. 

The  present  JAN  Standard  105  is  simply  the  Navy's  "Appendix  X" 
(General  Specifications  for  Inspection  of  Material,  Appendix  X,  Standard 
Sampling  Inspection  Tables  for  Inspection  by  Attributes,  issued  by  the  Navy 
Department,  April  1946)  with  a new  cover,  converting  it(|into  "JAN-STD  105, 
Sampling  Inspection  Tables  for  Inspection  by  Attributes)  and  bearing  such 
changes  in  the  general  instructions  as  were  necessary  to  convert  it  to  a 
joint  Army -Navy -Air  Force  document  pending  the  promulgation  of  a revised 
standard.  Various  changes  in  the  present  JAN  Standard  105  have  seemed 
necessary,  or  at  least  desirable,  in  the  light  of  the  needs  of  the  various 
inspection  agencies  of  the  Department  of  Defense,  and  their  experiences 
with  "Appendix  X",  the  present  standard,  and  related  documents.  Among  such 
changes  are: 

1)  Replacement  of  the  AQL  ranges  used  as  column  headings  of  the 
tables  by  single  AQL  values.  The  demand  for  this  change  was  especially 
strong  from  inspection  personnel  at  the  operating  level. 

2)  Provision  of  tables  ( i .e .,' acceptance -sampling  plans)  for 
larger  (l.e.  "looser")  AQL  values.  The  largest  AQL  accommodated  by"Appen- 
dix  X"  is  6.4$,  whereas  in  the  Quartermaster  Corps,  Signal  Corps,  and 
elsewhere,  inspection  by  attributes  may  be  called  for  in  cases  where  AQL's 
of  15$  or  20$  are  characteristic  of  the  product  concerned. 

3)  Inclusion  of  tables  for  inspection  on  a defects-per-unit 
basis  as  well  as  for  inspection  on  a percent  defective  (i.e.,  percent- 
having -at-least-one-defect)  basis.  Experience  of  the  Quartermaster  Corps, 
Signal  Corps,  Army  Ordnance  Inspection  services,  and  others,  with  defects- 
per-unit  inspection  has  seemed  to  indicate  that  the  indoctrination  and 
training  of  inspectors  at  the  operating  level  would  be  simplified  by  merging 
percent-defective  and  defects-per-unit  inspection  procedures  in  a single 
standard . 

4)  Replacement  of  Acceptable-Quality  Level  (AQL)  by  Lot-Quality 
Standard  (IQS)  . The  concept  of  an  Acceptable -Quality  Level  as  originally  in- 
troduced referred  to  "the  maximum  percent  defective  which  can  be  considered 
satisfactory  as  a process  average",  and  was  used  in  this  sense  in  standard 
inspection  procedures  issued  during  World  War  II  by  the  Army  Ordnance  De- 
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partment,  Quartermaster  Corps,  and  Signal  Corps.  It  thus  represented  an 
engineering  standard  of  permissible  non -perfection  adopted  in  recognition 
of  the  fact  that  to  demand  perfection  in  purchased  materials  was  unrealis- 
tic, or  at  least.  Impractical  and  uneconomic j and  was  related  to  long-run 
experience,  that  is,  to  the  process  average.  The  Navy's  "Appendix  X",  on 
the  other  hand,  introduced  the  same  term  and  'symbol  (AQL)  as  the  "percen- 
tage of  defective  items  in  an  inspection  lot  such  that  the  sampling  plan 
will  result  In  the  acceptance  of  95$  of  submitted  Inspection  lots  containing 
that  percentage  of  submitted  material".  The  AQL  so  defined  was  thus  a 
property  of  the  acceptance-sampling  procedure  concerned,  and  related  to  the 
acceptance  of  lots  by  the  procedure.  Much  confusion  has  resulted  from  the 
two-fold  difference  between  these  two  definitions  of  the  same  term  and  sym- 
bol. Administrative  and  operating-level  inspection  personnel  seem,  in  the 
main,  to  find  the  engineering-standard  approach  more  natural,  and  easier  to 
comprehend  and  employ.  Legal  counsel  advise  that  the  contractual  relation- 
ship between  the  purchaser  and  the  supplier  generally  relates  to  the  pur- 
chase of  lots,  and  that  the  quality  requirement  must  therefore  be  expressed 
in  terms  of  lots  not  process  averages.  Accordingly,  an  engineering  standard 
representing  the  contracted-for  quality  of  lots  seems  to  be  needed. 

5)  Provision  of  an  orderly  reduction  of  the  Producer's  Risk  with 
increasing  lot  size.  In  view  of  the  somewhat  different  practical  circum- 
stances attending  the  disposition  of  a large  rejected  lot  and  a small  re- 
jected lot,  the  risk  of  erroneously  rejecting  a lot  of  contracted-for  quality 
or  better  should,  It  seems,  be  less  for  large  than  for  small  lots.  Also, 
with  lot-by-lot  acceptance  sampling  It  is  generally  not  possible,  to  pro- 
vide the  same  protection  to  the  purchaser  against  acceptance  of  small  lots 
of  very  poor  quality  as  against  large  lots  of  such  quality.  Some,  often 
considerable.  Improvement  in  this  situation,  from  the  purchaser's  viewpoint, 
can  be  achieved  by  increasing  slightly  for  small  lots  the  producer's  risk 
of  having  lots  of  contracted-for  quality  erroneously  rejected  by  the  accep- 
tance-sampling procedure. 

As  the  quarter  drew  to  a close,  a final  draft  of  the  Sub-Committee' 
proposed  revision  of  JAN  Standard  105  neared  completion.  Thus  the  end  of 
this  the  first  leg  of  the  journey  toward  the  ultimate  goal  - a Government -wide 
standard  for  acceptance-sampling  by  attributes  - seems  to  be  in  sight. 


GLOSSARY  OF  STATISTICAL  ENGINEERING  TERMINOLOGY 
Project  11.3/2-48-3 

Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  INACTIVE. 


BIBLIOGRAPHY  AND  GUIDE  TO  STATISTICAL  LITERATURE 
Project  11.3/2-49-1 

Origin:  NBS 

Full  project  description  appears  in  Jan-Mar  1949  issue. 


Status:  INACTIVE. 


Status  of  Projects 


47 


GUIDE  TO  TABLES  OF  NORMAL  PROBABILITY  INTEGRAL 
Project  11.3/2-49-3 

Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 
Status:  INACTIVE. 


3 . Statistical  Services 


STATISTICAL  STUDIES  OF  CLINICAL  THERMOMETER  TESTING 
Project  11.3/31-47-3 

Origin:  NBS,  Section  3.1  Authorized  7/1/47 

Manager:  J.  M.  Cameron  Completed  9/30/49 

Objective:  To  determine  whether,  and  In  what  form,  acceptance 

sampling  is  feasible  in  inspection  and  testing  of  clinical  thermometers. 

Background:  Various  Government  agencies  (e.g..  Veterans  Administra- 

tion) that  purchase  clinical  thermometers  under  contracts  referencing 
Federal  Specification  GG-T-311;  Thermometers,  Clinical,  use  the  Bureau 
as  the  inspection  and  testing  agency  for  these  thermometers.  Although  the 
aforementioned  specification  includes  a provision  (Par.  F-3)  for  acceptance 
or  rejection  of  an  entire  delivery  on  the  basis  of  the  characteristics  and 
performance  of  a sample  of  thermometers  from  the  delivery,  many  contracts 
reference  the  paragraph  (F-2)  calling  for  complete  inspection  and  testing  of 
100$  of  the  thermometers  of  a delivery.  At  times  backlogs  of  serious  pro- 
portions have  resulted  from  this  practice.  The  present  project  is  intended 
to  explore  ways  and  means  of  reducing  the  volume  of  inspection  and  testing 
without  significant  loss  of  protection  to  the  purchasers. 

Comments:  Work  on  an  additional  report,  "Analysis  of  the  reproduci- 

bility of  the  tests  for  accuracy  and  consistency  of  readings,"  initiated 
under  this  project,  was  discontinued  and  will  be  renewed  under  another 
project,  project  11.3/31-49-1. 

Status:  COMPLETED.  Two  cooperative  experiments  carried  out  at  the 

Johns  Hopkins  Hospital  to  determine  (l)  the  reproducibility  of  clinical 
temperature  readings  under  customary  hospital  conditions  and  procedures 
and  (2)  the  correlation  of  the  NBS  mechanical  hard-shaker  tests  with  manual 
shaking  results  of  doctors  and  nurses,  were  analyzed  and  evaluated  and  the 
results  reported  by  letter  to  Dr.  Rider  of  the  Johns  Hopkins  Hospital. 

Other  work  bearing  on  this  project  was  transferred  to  project  11.3/31-49-1 
under  which  more  suitable  data  recently  became  available. 
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WOOL  CONTENT  OF  BLANKETS 
Project  11.5/J1-47-6 

Origin:  NBS  and  Division  of  Statistical  Standards,  Bureau  of  the  Budget. 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Work  on  this  project  from  March  to  September 

1949,  took  the  form  of  active  participation  by  Dr.  Eisenhart  in  the  work 
of  a Task  Group  on  Revision  of  Army  Material  and  Test  Method  Specifica- 
tions for  Textiles,  set  up  under  Sub-Committee  B-5  on  Presentation  and 
Interpretation  of  Data,  of  Committee  D-13  on  Textiles,  of  the  American 
Society  for  Testing  Materials,  in  response  to  a request  from  the  Sub- 
committee on  Fibers  and  Fabrics  of  the  National  Research  Council's  Ad- 
visory Board  on  Quartermaster  Research  and  Development. 

The  membership  of  the  Task  Group  consists  of  Mr . 0.  P.  Beckwith 
(Alexander  Smith  and  Sons  Carpet  Company),  chairman,  Stanley  Backer  (Re- 
search and  Development  Branch.  Quartermaster  Corps) , Churchill  Eizenhart 
(National  Bureau  of  Standards),  John  C.  Hintermaier  (Cluett  Peabody  and 
Company),  Arthur  Lipman  (Philadelphia  Quartermaster  Depot),  David  H. Schwartz 
(Quartermaster  Corps  Inspection  Service),  and  Gerald  Winston  (Research 
and  Development  Brand).,  Quartermaster  Corps)  . Each  member  serves  solely 
as  an  expert  in  his  field,  not  as  an  official  representative  of  his  or- 
ganization . 

During  recent  months  there  had  been  numerous  instances  of  dis- 
agreement over  the  interpretation  of  physical-  and  chemical-test  require- 
ments of  Quartermaster  Corps  textile-materials  and  textile -test-methods 
specifications.  Study  of  the  situation  revealed  that  there  was  a dif- 
ference between  the  intent  of  the  Research  and  Development  specification 
writer  and  the  interpretation  of  the  Inspection  Services  and  the  contrac- 
tors. The  Task  Group  was  charged  with  the  clarification  of  such  spec- 
ifications so  as  to  leave  no  room  for  misinterpretation,  and  to  eliminate 
any  inconsistencies  between  test-methods  specifications  and  material  spec- 
ifications. In  this  connection,  the  following  ambiguities  or  omissions 
were  found  to  be  characteristic  of  existing  specifications  of  the  types 
under  consideration:  1)  ambiguity  on  whether  the  physical  and  chemical 

requirements  applied  to  the  average  of  all  the  units  in  a lot  or  to  each 
unit  individually,  2)  omission  of  information  on  the  quantity  of  cloth 
constituting  an  individual  "unit"  of  production  for  test  purposes, 

5)  ambiguity  regarding  the  number  of  measurements  of  a given  property  or 
characteristic  to  be  made  on  each  unit,  4)  omission  of  information  on  the 
consistency  of  repeated  determinations  of  a particular  property  made  in 
accordance  with  a prescribed  method  under  controlled  conditions  by  a single 
operator-equipment  combination  and  by  different  operator-equipment  combin- 
ations, and  ambiguity  in  the  relation  of  such  information,  in  conjunction 
with  the  variability  of  the  material  itself,  to  the  required  number  of  de- 
terminations per  unit,  and  5)  omission  of  an  engineering  standard  of 
permissible  nonperfection,  representing  contracted-for  lot  quality.  In 
recognition  of  the  fact  that  to  demand  perfection  (i.e.  no  defective  units, 
no  defects  per  unit)  is  unrealistic,  or  at  least  impractical  and  uneconomic. 

A proposed  revision  of  the  Quartermaster  Corps  textile-test- 
method  specification  and  a model  material  specification  embodying  the  above 
considerations  were  prepared  by  the  Task  Group  and  submitted  late  in 
September  1949  to  the  Sub-Committee  on  Fibers  and  Fabrics  for  its  con- 
sideration . 
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EFFECT  OF  GASOLINE  AND  OIL  ADDITIVES  ON  CARBON 
AND  GUM  FORMATION 
Project  11.3/31-47-8 


Origin:  NBS,  Division  3 Authorized  7/1/47 

Managers:  J.  M.  Cameron,  W.  J.  Youden  Completed  9/30/49 

Objective:  To  advise  and  assist  personnel  of  Section  3-5,  NBS, 
with  the  statistical  aspects  of  the  planning  and  conduct  of  experiments 
to  determine  the  effect  of  gasoline  and  oil  additives  on  carbon  and  gum 
formation  of  engines. 

Background:  By  utilization  of  recent  advances  in  the  techniques 

of  statistical  inference  and  the  principles  of  experimental  design,  it 
is  expected  that  economy  and  increased  efficiency  will  be  effected  in  this 
research  and  testing  program.  The  experiment  involves  the  testing  of 
over  20  additives  in  combination  with  a "control"  gasoline  and  oil  mix- 
ture on  80  similar  engines,  the  problem  being  to  design  the  most  efficient 
experiment  to  determine  the  performance  of  the  various  additives. 

Status:  COMPLETED.  An  interoffice  memorandum  summarizing  the 

statistical  aspects  of  the  tests  was  sent  to  Division  3- 


STATISTICAL  STUDY  OF  THE  FLOW  OF  CERTAIN  STOCKROOM  ITEMS 

Project  11.3/31-48-4 


Origin:  NBS,  Section  M.4 

Full  project  description  appears  in  Jan-Mar  1949  issue. 
Status:  INACTIVE. 


STATISTICAL  ANALYSIS  OF  THERMOMETRIC  MEASUREMENTS 
Project  11.3/31-49-4 


Origin:  NBS  , Section  3-1 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Analysis  of  the  data  from  all  6 of  the  tests 

carried  out  by  the  Temperature  Measurements  Section  (3-1)  was  completed, 
and  the  results  were  presented  informally  to  that  section. 


WOMEN'S  BODY  MEASUREMENT  STUDY 
Project  11.3/31-49-5 

Origin:  NBS,  Section  12.2 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Certain  age  groupings  were  decided  upon  after 

a study  -of  the  frequency  distributions  by  single  year  of  age  and  for  seven 
control  dimensions  (stature,  hip  girth,  bust  girth,  waist  girth,  abdominal 
extension  girth,  abdominal  extension  arc,  neck  to  bust  length) . Stature- 
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hip  girth  bivariate  frequency  distributions  were  then  obtained  for  each 
of  the  age  groupings.  From  these,  recommendations  were  made  for  four 
stature  groups  (short-short,  short,  regular,  tall). 

During  the  quarter  several  conferences  of  the  Washington  area 
group  were  held  and  one  conference  was  held  in  New  York  City  with  re- 
presentatives of  the  Mail  Order  Association  of  America,  to  whom  recom- 
mendations were  presented  and  plans  developed  for  continuation  of  the 
project . 


EXPERIMENTAL  ARRANGEMENTS 
Pro  ject 


FOR  MEASUREMENT  OF  RADIOACTIVITY 
11.3/31-50-1 


Origin:  NBS,  Section  4.4 

Manager:  W.  J.  Youden 


Authorized  8/1/49 
Completed  9/30/49 


Objective:  To  review  the  statistical  aspects  of  the  problem  of 

comparing  standards  of  radioactivity. 


Background:  Mr.  H.  H.  Seliger  pointed  out  that  extreme  care  was 

necessary  to  avoid  drifts  and  changing  conditions  when  using  Geiger 
counters  for  comparing  the  radioactivity  of  different  sources.  Cus- 
tomary practice  achieves  this  by  making  every  second  observation  a de- 
termination with  the  primary  standard. 

Status:  COMPLETED.  Experimental  trials  showed  that  a Latin 

Square  design  accomplished  the  same  purpose  of  controlling  the  effects 
of  changing  conditions  with  a marked  reduction  in  the  number  of  deter- 
minations on  the  standard. 


PRECISION  MEASUREMENTS  ON  STANDARD  TEMPERATURE  SOURCES 
Project  11.3/31-50-2 

Origin:  NBS,  Section  5.0  Authorized  8/1/49 

Manager:  W.  J.  Youden 

Objective:  To  ascertain  the  precision  of  measurements  directed 

to  the  determination  of  the  reproducibility  of  temperatures  set  by  the 
freezing  point  of  chemical  substances . 

Background:  A previous  study  by  F.  W.  Schwab  and  Edward  Wichers 

(NBS  RP  1647)  on  the  freezing  temperature  of  benzoic  acid  is  being 
extended  to  other  chemicals.  The  reproducibility  of  these  temperature 
sources  is  apparently  close  to  the  limits  detectable  by  the  most  re- 
fined methods  of  temperature  measurement.  It  is  necessary  to  devise 
methods  to  determine  whether  the  small  differences  found  between  re- 
peated measurements  may  be  ascribed  to  the  temperature  source  or  the 
method  of  measuring  the  temperature . 

Status:  NEW.  A schedule  for  the  intercomparison  of  sources  and 

thermometers  was  developed  which  increased  the  precision  of  the  com- 
parison of  sources  and  thermometers.  Consideration  is  being  given  to 
means  of  determining  whether  the  source  of  the  temperature  or  the  method 
of  measuring  the  temperature  is  the  limiting  factor  in  the  observed 
result . 
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PROBABILITY  THEORY  OF  A CERTAIN  DOUBLE  BRANCHING  SYSTEM 
Project  11.3/51-50-5 


Origin:  NBS,  Section  4.5  Authorized  9/1/49 

Manager:  J.  Lietlein 

Objective:  To  study  the  probability  distribution  of  two  types  of 

particles  in  a certain  multiplicative  system,  with  particular  reference 
to  developing  statistical  tests  of  the  "crossover"  probability  based 
on  observed  frequencies  . 

Background:  Each  of  an  initial  set  of  similar  particles  (type  A) 

splits  into  two,  and  each  of  the  two  new  particles  has  probability  P (the 
"crossover"  probability  of  the  system)  of  being  a new  type,  B,  and 
probability  (l-P)  of  remaining  of  type  A.  If  the  new  particle  is  of 
type  A,  it,  in  turn,  can  produce  both  types  A and  B,  according  to  the  same 
probability  law,  but  if  of  type  B it  can  produce  only  type  B.  The  suc- 
cessive generations  of  A and  B particles  thus  produced  constitute  an  im- 
portant type  of  multiplicative  (branching)  process  in  two  variables  which 
does  not  appear  to  have  been  specifically  treated  heretofore.  This  prob- 
lem has  been  communicated  by  Dr.  Fano  in  connection  with  biophysical  in- 
vestigations . 

Status:  NEW.  About  60^  completed.  Explicit  values  for  several 

of  the  moments  of  the  distribution  of  both  types  of  particles  in  the  nth 
generation  have  been  evaluated  and  other  features  of  the  distribution 
were  in  process  of  investigation.  Everett  and  Ulam  have  developed  a 
theory  which  implicitly  contains  these  results,  but  which,  because  of 
its  completely  general  formulation,  is  in  terms  not  practicable  to  apply 
in  the  present  problem.  It  therefore  was  deemed  advisable  to  attack  the 
problem  directly  by  adapting  and  developing  some  of  the  theory  outlined 
in  a paper  by  T.  E.  Harris  dealing  with  branching  systems  of  only  one 
type  of  particle,  thereby  leading  more  directly  and  efficiently  to  the 
desired  results. 


STATISTICAL  SERVICES  TO  ARMY  FIELD  FORCES  TEST  BOARDS 
Project  11.3/32-49-1 

Origin:  Logistics  Division,  General  Staff,  U.S.A  Authorized  7/1/48 

Sponsor:  Research  and  Development  Group,  General  Staff, 

U . S . Army 

Manager:  W.  J.  Youden 

Objective:  To  provide  consultation  services  to  the  Army  Field 

Forces  Test  Boards  and  on  other  projects  of  interest  to  the  Logistics 
Division . 

Background:  The  Test  Boards  are  concerned  with  planning  experimen- 

tal field  tests  of  Army  equipment  and  weapons  and  the  evaluation  of  data 
where  the  nature  of  the  tests  limits  the  number  of  the  objects  tested 
to  small  numbers.  It  is  desirable  to  make  use  of  modern  statistical  tech- 
niques in  the  conduct  of  the  test  programs. 

Comments:  This  is  a continuing  project.  Visits  are  made  to  the 

Boards  either  on  request  or  to  make  available  the  results  of  studies 
originating  as  a consequence  of  previous  visits.  On  occasion  extended 
visits  are  made  (in  one  case  nearly  a month)  for  the  purpose  of  becoming 
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familiar  with  the  problems  of  the  Test  Boards  by  working  alongside  of  the 
officer  personnel  on  current  projects. 

Status:  CONTINUED.  (l)  Criteria  and  brief  tables  were  provided  for 

comparing  the  aiming  errors  of  one  system  with  those  of  another  system, 
and  the  aiming  error  of  a given  system  with  a prescribed  performance  when 
only  a few  measurements  are  in  hand.  (2)  In  the  event  that  one  of  the 
aiming  errors  appears  unusually  large  and  is  therefore  possibly  a gross 
error  it  was  pointed  out  the  median  error  has  certain  advantages  or  alter- 
natively that  a statistical  criterion  is  available  to  assist  in  deter- 
mining whether  the  measurement  in  question  may  be  considered  a gross  error. 
(3)  These  methods  were  extended  to  errors  measured  in  space  as  well  as  In 
two  dimensions.  (4)  The  investigation  of  aiming  errors  over  a range  of 
conditions  when  only  very  limited  numbers  of  firings  are  made  presents  a 
particularly  difficult  problem.  A technique  was  proposed  for  this  situation. 


SOLID  PROPELLANT  DEVELOPMENT  TESTS 
Project  11.3/32-49-2 

Origin:  ORDTU,  Research  and  Development  Division,  Ordnance  Dept., 

U . S . Army 

Sponsor:  Research  and  Development  Group,  General  Staff,  U.  S.  Army 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  INACTIVE 


PROBABILITY  STUDIES  OF  MISSILE  EFFECTIVENESS 
Project  11.3/33-49-1 

Origin:  Bureau  of  Ordnance,  Navy  Department,  via  Electronics  Division, 

NBS 

Sponsor:  Bureau  of  Ordnance,  Department  of  the  Navy 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  With  the  aid  of  Professor  A.  H.  Bowker, 

acting  chairman  of  the  Department  of  Mathematical  Statistics,  Stanford 
University,  who  served  as  a special  consultant  for  three  weeks,  the 
basic  assumptions  and  probability  models  were  set  up.  The  various  data 
necessary  for  the  solution  of  the  problem  were  listed,  and  action  was 
taken  on  their  procurement . 


TECHNICAL  SYMBOLISM  AND  TERMINOLOGY 
Project  11.3/4-50-1 

Origin:  NBS  Educational  Committee  Authorized  9/1/49 

Manager:  Lola  S.  Deming 

Objective:  To  offer  an  in-hours  "editorial  course"  designed  es- 

pecially for  typists  and  clerks  who  are  responsible  for  transcribing, 
typing,  proofing,  and  editing  scientific  material.  Topics  covered  to 
Include:  Greek  alphabet  (recognition,  pronunciation,  writing,  and  typing); 
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graphs,  charts,  and  tables  (scales,  legends,  general  format,  etc.); 
exponents,  subscripts,  radicals,  other  mathematical  and  chemical  sym- 
bols; vocabulary  and  spelling  of  common  scientific  terms  and  proper 
names;  rounding  (rules  and  practices);  significant  figures. 

Background:  Because  the  preparation  of  a large  proportion  of  man- 

uscripts at  the  Bureau  involves  transcribing,  typing,  and  editing  of 
mathematical  and  scientific  material,  a course  of  this  nature  seemed 
useful.  It  was  offered  for  the  first  time  last  year  and  was  received 
with  sufficient  enthusiasm  to  cause  the  Educational  Committee  to  ask 
to  have  it  repeated  this  year  and  extended  to  two  terms  to  Include  more 
material  with  several  sessions  conducted  by  a chemist  to  cover  chemical 
terminology  and  vocabulary. 

Comments:  This  course  is  listed,  with  the  above  title,  and  de- 

signated Flll.l,  on  page  17  of  the  Announcement  of  Courses  for  19^9-1950 
of  the  Graduate  School  of  the  National  Bureau  of  Standards. 

Status:  NEW.  In  the  absence  of  appropriate  published  text  material, 

dittoed  notes  were  made  for  distribution  to  the  class  and  to  serve  later 
as  the  backbone  of  a manual  for  the  guidance  of  NBS  clerks  and  typists 
in  preparing,  proofing,  and  editing  technical  manuscripts. 


IV.  Machine  Development  Laboratory 
(Section  11.4) 


1.  Design  and  Construction  of  Automatic  Digital  Computing  Machines 


The  machine  design  and  construction  projects  are  being  performed 
in  cooperation  with  the  Electronics  Division  of  the  Bureau. 


THE  BUREAU  OP  THE  CENSUS  COMPUTING  MACHINE 
Project  11.4/21-47-1 

Origin:  The  Bureau  of  the  Census 

Sponsor:  The  Bureau  of  the  Census 

Full  project  description  appears  in  Apr-Jun  1949  Issue. 

Status:  CONTINUED.  Two  check  points  of  the  UNIVAC  contract  were 

successfully  passed  by  the  contractor,  one  relating  to  engineering 
design  (the  BINAC  check  point)  and  the  other  concerned  with  auxiliary 
equipment  (the  UNISERVO  check  point).  UNIVAC  design  modifications, 
recommended  by  the  contractor  to  expedite  construction,  were  found  to 
decrease  the  UNIVAC  operational  speed  by  15  to  20  per  cent,  but  were 
considered  acceptable . 


THE  NAVY  COMPUTING  MACHINE 
Project  11.4/22-47-2 

Origin:  Mathematics  Branch,  Office  of  Naval  Research 

Sponsor:  Office  of  Naval  Research 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  logical  design  of  this  computer  is  essen- 

tially complete.  Instruction  codes  have  been  made  available  by  the 
Raytheon  Company  for  analysis  by  the  programing  group  of  this  laboratory 
Mercury-line  storage  and  tape -handling  units  were  placed  under  construc- 
tion during  the  quarter,  and  in  addition  fairly  extensive  prototyping 
of  computer  components  was  under  way. 


AIR  MATERIEL  COMMAND  COMPUTING  MACHINE 
Project  11.4/23-49-1 

Origin:  Air  Materiel  Command.  USAF 

Sponsor:  Office  of  Air  Research,  AMC , USAF 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Computer  proposals  by  the  Engineering  Research 

Associates  were  proposals  of  machines  whose  cost  would  exceed  the  funds 
available  for  the  project.  General  Electric  cost  estimates  were  lower 
but  were  estimates  made  by  design  engineers  and  had  not  been  submitted 
to  the  General  Electric  Production  Department.  Plans  for  advertising 
for  bids  for  this  computer  had  been  made  by  the  end  of  the  quarter. 
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THE  AIR  COMPTROLLER'S  COMPUTING  MACHINE 

Project  11.4/24-47-3 

Origin:  Office  of  the  Air  Comptroller,  USAE 

Sponsor:  " " " " 

Pull  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  This  computer  is  being  constructed  by  the 

Eckert-Mauchly  Computer  Corporation.  The  contract  is  in  the  form  of  a 
supplement  to  the  Census  UNIVAC  contract.  The  status  Is  the  same  as  given 
for  project  11.4/21-47-1. 


THE  NBS  INTERIM  COMPUTER 
Project  11.4/24-49-1 

Origin:  KBS 

Sponsor:  Air  Comptroller's  Office,  USAP 

Pull  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  logical  design  has  been  frozen  . The 

memory  unit  has  been  fully  engineered  and  sufficient  prototypes  of 
Important  components  have  been  built  and  tested  to  provide  the  contrac- 
tor with  the  data  necessary  for  him  to  proceed  with  construction. 
Construction  of  components  has  been  under  way  by  the  Electronic  Computer 
Section  of  the  Bureau.  It  is  expected  that  the  computer  will  be  as- 
sembled and  under  test  early  in  1950. 


ARMY  MAP  SERVICE  COMPUTING  MACHINE 
Project  11.4/25-49-1 

Origin:  Army  Map  Service,  U.S.A. 

Sponsor:  " " " 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  Army  Map  Service  Computer  construction 

contract,  with  the  Eckert-Mauchly  Computer  Corporation,  is  in  the  form 
of  a supplement  to  the  Census  UNIVAC  contract,  providing  for  delivery 
of  a duplicate  of  the  Census  UNIVAC  system.  The  status  is  the  same  as 
given  for  project  11.4/21-47-1. 


DEPARTMENT  OF  THE  ARMY  COMPUTER  DESIGN 
Project  11.4/26-49-1 

Origin:  Department  of  the  Army 

Sponsor:  " " 

Full  project  description  appears  In  Apr-Jun  1949  issue. 
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Status:  CONTINUED.  The  memory  units  are  under  construction  by 

the  contractor.  Little  direct  work  on  the  machine  was  done  during  the 
quarter  under  review;  however,  much  of  the  work  on  the  NBS  Interim 
Computer  was  applicable  to  this  project  (see  project  11.4/24-49-1) 


2 . Programing  Studies 


PROGRAMING  OF  PROBLEM^  FOR  SOLUTION  ON  AUTOMATIC  DIGITAL 

COMPUTING  MACHINES 
Project  11.4/3-47-4 

Origin:  Bureau  of  the  Census,  Department  of  the  Navy,  Department  of 

the  Air  Force,  and  Department  of  the  Army. 

Sponsors:  Bureau  of  the  Census,  Department  of  the  Navy,  Department 

of  the  Air  Force,  and  Department  of  the  Army. 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Emphasis  was  placed  upon  investigation  of 

the  utility  of  the  UNIVAC  as  a sorting  device.  A working  group,  con- 
sisting of  representatives  of  the  Bureau,  the  Bureau  of  the  Census, 
the  USAF  and  the  Navy,  was  organized  to  study  techniques  of  electronic 
sorting  and  to  evaluate  existing  and  proposed  devices  for  high-speed 
sorting . 


CODING  ON  THE  E.R.A.  COMPUTER 
Project  11.4/3-49-1 


Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  Programing  analyses  were  performed  in  con- 

nection with  the  evaluation  of  the  E.R.A.  proposals  for  an  Air  Materiel 
Command  computer  (see  project  11.4/23-49-1). 


CODING  RELATED  TO  THE  UNIVAC  SYSTEM 
Project  11.4/3-49-2 

Origin:  The  Bureau  of  the  Census 

Sponsor:  " " " 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  An  analysis  of  the  latest  Instruction  code 

submitted  by  the  Eckert-Mauchly  Computer  Corporation  for  the  UNIVAC  was 
made  in  order  to  determine  the  effect  of  the  suggested  changes  on  the 
programing  of  problems  for  the  computer  and  the  relative  efficiency  of 
the  new  instructions. 
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CODING  RELATED  TO  THE  RAYTHEON  COMPUTER 
Project  11.4/3-49-5 

Origin:  Mathematics  Branch,  Office  of  Naval  Research 

Sponsor:  Office  of  Naval  Research,  USN 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  Raytheon  Manufacturing  Company  revised  the 

codes  and  design  of  Its  machine.  Some  programing  was  done  for  the  purpose 
of  comparing  the  speed  of  its  performance  with  that  of  the  UNIVAC . 


PERFORMANCE  TEST  FOR  THE  BINAC 
Project  11.4/3-49-4 

Origin:  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  test  routines  for  the  BINAC  which  had  been 

prepared  at  the  Bureau  were  rejected  by  the  Eckert-Mauchly  Corporation 
as  too  time-consuming,  and  they  submitted  two  of  their  own  tests  instead. 
These  tests  were  analyzed  and  that  one  chosen  which  more  closely  tested 
the  performance  specifications  desired  by  the  Bureau.  The  test  was  ap- 
plied, and  with  the  assistance  of  the  Census  Bureau  the  results  were 
checked  and  evaluated. 


3 . Technical  Reports  on  Computing  Machinery 


THE  MTAC  SECTION 
Project  11.4/4-47—1 

Origin:  Committee  on  High-Speed  Computing  of  the  National  Research 

Council 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  The  material  for  the  Automatic  Computing 

Machinery  section  of  the  January  1950  issue  of  Mathematical  Tables  and 
Other  Aids  to  Computation  was  assembled  and  edited.  A technical  paper 
entitled  "Magnetic  Drum  Storage  for  Digital  Information  Processing  Sys- 
tems," by  Arnold  A.  Cohen,  describes  a machine  constructed  by  Engin- 
eering Research  Associates,  Inc.  In  a discussion  article,  "Solution 
of  Differential  Equations  by  Recurrence  Relations,"  by  John  Todd,  the 
author  shows,  by  consideration  of  a very  simple  example,  that  recursive 
methods  for  the  solution  of  ordinary  differential  equations,  while 
very  attractive  from  the  point  of  view  of  mechanization,  must  be  used 
with  care.  The  example  is  discussed  in  some  detail  to  exemplify  a few 
of  the  concepts  used  in  modern  numerical  analysis,  such  as  truncation 
error,  rounding-off  error,  strict  estimates,  probabilistic  estimates, 
and  stability.  Attention  Is  drawn  to  work  of  L.  Collatz,  who  discussed 
similar  problems  for  partial  differential  equations,  and  to  recent 
work  of  L.  Fox  and  E.  T.  Goodwin,  who  have  developed  an  iterative 
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process  for  the  use  of  recursion  formulae  and  applied  it  successfully 
in  several  cases  . 


BIBLIOGRAPHY  ON  HIGH-SPEED  AUTOMATIC  COMPUTING  MACHINERY 

Project  11.4/42-49-2 

Origin:  Section  11.4,  NBS 

Full  project  description  appears  in  Apr-Jun  1949  issue. 

Status:  CONTINUED.  About  85$  completed. 


Lectures  and  Symposia 


HARTREE,  D. 

LANCZOS,  C 
HELMER,  0. 

ISAACS,  R. 
DONSKER,  M. 

GRUEN,  H. 


BECKENBACH 
OSTROVSKI , 

ROSSER,  J. 
SZASZ , 0. 
TODD,  J. 


HUSKEY,  H. 
LANCZOS , C 


Numerical  Analysis  Colloquium  Series 
(Los  Angeles,  California) 


R.  (The  Cavendish  Laboratory,  Cambridge,  England).  (l)  "Forced 
oscillation  of  non-linear  systems."  July  6,  1949.  (2)  "Some 

electron  space-charge  calculations."  July  f,  194-9.  (3)  "A 

method  for  the  numerical  integration  of  first-order  differen- 
tial equations."  July  8,  1949. 

"The  calculation  of  eigenvalues."  July  25  ■,  1949- 

(Rand  Corporation).  "Methods  in  the  theory  of  games."  August 

1,  1949. 

(Rand  Corporation).  "Monodiffric  functions ."  August  8,  1949- 

"Calculation  of  eigenvalues  of  the  Schrodinger  equation  by 
Monte  Carlo  methods."  August  15,  1949. 

(U.C.L.A.  and  National  Bureau  of  Standards.)  "Properties  of 
light  nuclei."  August  22,  1949. 


Papers  and  Invited  Talks 
Presented  by  Members  of  the  Staff 
at  Meetings  of  Outside  Organizations 


1)  Papers  presented  at  the  meeting  of  the  American  Mathematical 
Society  held  at  the  University  of  Colorado,  August  30- 
September  2,  1949: 

E.  F.  "A  class  of  mean  value  functions." 

A.  M.  "On  Theordorsen ' s and  Garrick's  method  of  conformal 

mapping . " 

B.  "Real  roots  of  Dirichlet  L-series." 

"On  the  Gibbs'  phenomenon  for  Euler  means." 

'The  condition  of  certain  matrices."  Presented  by  title. 

2)  Papers  presented  at  A Symposium  on  Large-Scale  Digital  Com- 
puting Machinery  held  at  Harvard  University  under  the 
sponsorship  of  the  Navy  Department  Bureau  of  Ordnance 

and  Harvard  University,  September  13-16,  1949: 

D.  "Semi-automatic  Instruction  on  the  Zephyr." 

"An  iteration  method  for  the  solution  of  the  eigenvalue 
problem  of  linear  differential  and  integral  operators." 
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MILNE,  W.  E.  "Numerical  methods  associated  with  Laplace's  equation." 

(3)  Invited  addresses  presented  at  meetings  of  the  Peripatetic 
Society. 

HUSKEY,  H.  D.  "The  National  Bureau  of  Standards  electronic  digital 

computer  program."  At  the  University  of  Southern  Cal- 
ifornia, August  22,  1949. 

KAC,  M.  "Zeros  of  random  functions."  At  the  University  of  California 
at  Los  Angeles,  July  25,  1949. 

WARSCHAWSKI , S.  E.  "On  the  behavior  of  the  mapping  function  in  conformal 
mapping."  At  the  University  of  Southern  California, 
August  22,  1949. 

(4)  Invited  talks  given  at  the  Rand  Corporation,  Santa  Monica, 
California . 

HESTENES,  M.  R.  "Problem  of  Bolza  with  applications  to  paths  of  least 
time."  September  13,  1949. 

KAC,  M.  "Deviations  between  theoretical  and  empirical  distributions." 
August  2,  1949. 

(5)  Other  invited  talks: 

MILNE,  W.  E.  "A  new  method  of  numerical  solution  of  differential  equa- 
tions." Presented  at  the  National  Bureau  of  Standards, 
Washington,  D.  C.,  September  9,  1949- 

0STR0WSKI,  A.  M.  "On  the  continuity  of  roots  of  algebraic  equations." 

Presented  to  the  Mathematics  Department  of  the  University 
of  Maryland,  September  22,  1949. 


Publication  Activities 


1 . Publications  which  appeared  during  the  quarter 

1.1  Mathematical  Tables 

(None) 

1.2  Manuals,  Bibliographies  and  Indices 

( None) 

1. J  Technical  Papers 

(1)  Formulas  for  the  percentage  points  of  the  distribution  of  the 
arithmetic  mean  in  random  samples  from  certain  symmetrical  uni- 
verses. Uttam  Chand . NBS  J.  Res.  43_,  No.  1,  79-80  (July  1949). 
Available  as  RP2007  from  the  Superintendent  of  Documents,  Gov- 
ernment Printing  Office,  Washington  25,  D.C.  at  10  cents. 

(2)  Probability  center  lines  for  standard  deviation  and  range 
charts.  C.  Eisenhart.  Industrial  Quality  Control,  VI,  No.  1, 
24-26  (July  1949).  Reprints  available. 

2.  Manuscripts  in  the  Process  of  Publication  September  30,  1949 

2.1  Mathematical  Tables 

(1)  Tables  of  the  binomial  probability  distribution.  NBS  Applied 
Mathematics  Series  6.  Being  printed  by  the  Government  Printing 
Office  . 

(2)  Tables  to  facilitate  sequential  tests.  NBS  Applied  Mathematics 
Series  7.  Being  printed  by  the  Government  Printing  Office. 

(3)  Tables  of  powers  of  complex  numbers.  NBS  Applied  Mathematics 
Series  8 . Being  printed  by  the  Government  Printing  Office . 

Exact  values  of  zn  for  z = x + iy,  x,  y = 0(l)l0,  n = 0(l)25; 
exact  values  of  xn  for  x = 0(1)10,  n = 0(1)25. 

(4)  Tables  of  the  Chebyshev  polynomials  Sn(x)  and  Cn(x)  . NBS  Ap- 
plied Mathematics  Series  9*  Being  printed  by  the  Government 
Printing  Office  . 

(5)  Tables  of  Bessel  functions  YQ(z)  and  Yp(z)  for  complex  arguments. 
Being  printed  by  Gblumbia  University  Press. 

(6)  Tables  relating  to  the  Mathieu  functions.  Being  printed  by  the 
Columbia  University  Press. 

(7)  Table  of  coefficients  in  numerical  integration  formulae.  NBS 
Mathematical  Table  22.  To  be  reissued.  Being  printed  by  the 
Government  Printing  Office.  (Originally  published  in  J.  Math. 
Phys.  XXII,  No.  2,  49-50,  June  1943.) 

(8)  Table  of  Fourier  coefficients.  NBS  Mathematical  Table  23.  To 
be  reissued.  Being  printed  by  the  Government  Printing  Office. 
(Originally  published  in  J.  Math.  Phys.  XXII,  No.  3,  136-147, 
Sept.  1943.) 
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2.J  Technical  Papers 


(1) 

Asymptotic  expansions  of  spheroidal  wave  functions.  M .Abramowitz . 
Accepted  for  publication  in  the  Journal  of  Mathematics  and  Physics 

(2) 

Tables  of  integrals  of  Struve  functions.  M.  Abrainowitz.  Accepted 
for  publication  in  the  Journal  of  Mathematics  and  Physics. 

(3) 

A class  of  mean  value  functions.  E.  F.  Beckenbach.  Submitted  to 
the  American  Mathematical  Monthly. 

(4) 

Metric  differential  geometry.  E.  F.  Beckenbach.  Submitted  to  the 
Mathematics  Magazine. 

(5) 

On  subordination  in  complex  theory.  E.  F.  Beckenbach  and 
E.  W.  Graham.  Submitted  to  the  Bulletin  of  the  American  Mathe- 
matical Society;  also  to  appear  in  the  "Proceedings  of  a sym- 
posium on  the  construction  and  applications  of  conformal  mapping," 
to  be  published  by  the  National  Bureau  of  Standards  . 

(6) 

Recurrent  determinants  of  orthogonal  polynomials.  Part  I:  Legen- 
dre and  ultraspherical  polynomials.  E.  F.  Beckenbach,  W.  Seidel 
and  0.  Szasz.  Submitted  to  the  Duke  Mathematical  Journal. 

(7) 

Table  of  modified  Bernoulli  polynomials.  G.  Blanch  and  R.  Siegel. 
Accepted  for  publication  in  the  NBS  Journal  of  Research. 

(8) 

A "Simpson's  rule"  for  the  numerical  evaluation  of  Wiener's 
integrals  in  function  space.  R.  H.  Cameron.  Submitted  to  the 
Duke  Mathematical  Journal. 

(9) 

The  application  of  statistical  procedures  to  the  preparation  of 
industrial  specifications  and  acceptance  procedures.  J.  H.  Curtis 
To  appear  in  the  Proceedings  of  the  International  Statistical  Con- 
ferences . 

(10) 

Acceptance  sampling  by  variables,  with  special  reference  to  the 
case  in  which  quality  is  measured  by  average  or  dispersion. 

J.  H.  Curtiss.  To  appear  in  a special  supplement  to  the  Journal 
of  the  American  Statistical  Association. 

(11) 

A sampling  method  for  determining  the  lowest  eigenvalue  and  the 
principal  eigenfunction  of  Schrddinger ' s equation.  M.  D.  Donsker 
and  M.  Kac . Accepted  for  publication  in  the  NBS  Journal  of  Re- 
search . 

(12) 

Exact  particle  trajectories  for  nonviscous  flow  in  a plane  with 
a constant  coriolis  parameter.  G.  E.  Forsythe.  Accepted  for 
publication  in  the  Journal  of  Meteorology.  (Originally  submit- 
ted under  the  title.  Some  exact  solutions  of  the  equations  of 
non-viscous  flow  in  a plane.) 

(13) 

Solution  of  the  telegrapher's  equation  with  boundary  conditions 
on  only  one  characteristic.  G.  E.  Forsythe.  Accepted  for  pub- 
lication in  the  NBS  Journal  of  Research. 
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(14)  Numerical  integration  for  linear  sums  of  exponential  functions. 
R.  E,.  Greenwood.  Submitted  to  the  Quarterly  of  Applied  Mathe- 
matics . 

(15)  1949  Ephemeris  of  Jupiter's  ninth  satellite.  S.  Herrick.  To 
be  published  as  a leaflet  by  the  Astronomical  Society  of  the 
Pacific  . 

(16)  Some  elementary  problems  in  the  calculus  of  variations.' 

M.  R.  Hestenes.  Submitted  to  the  Mathematics  Magazine;  also 
to  be  published  as  a chapter  in  a survey  of  mathematics  edited 
by  Professor  Glenn  James  of  U.C.L.A. 

(17)  Quadratic  forms  in  Hilbert  space,  with  applications  in  the 
calculus  of  variations.  M.  R.  Hestenes.  Accepted  for  pub- 
lication in  the  American  Journal  of  Mathematics. 

(18)  Systems  of  extremals  for  the  simplest  isoperimetric  problem. 

M.  Karlin.  Submitted  to  the  Bulletin  of  the  American  Mathe- 
matical Society. 

(19)  The  remainder  in  linear  methods  of  approximation.  W.  E.  Milne. 
Accepted  for  publication  in  the  NBS  Journal  of  Research. 

(20)  A note  on  the  numerical  integration  of  differential  equations. 
W.  E.  Milne.  Accepted  for  publication  in  the  NBS  Journal  of 
Research . 

(21)  Numerical  determination  of  characteristic  numbers.  W.E. Milne. 
Accepted  for  publication  in  the  NBS  Journal  of  Research. 

(22)  Generalization  of  a theorem  of  Osgood  to  the  case  of  continuous 
approximation.  A.  M.  Ostrowski.  Submitted  to  the  Bulletin 

of  the  American  Mathematical  Society. 

(23)  On  some  generalizations  of  the  Cauchy -Frullani  integral. 

A.  M.  Ostrowski.  Accepted  for  publication  in  the  Proceedings 
of  the  National  Academy  of  Sciences. 

(24)  Note  on  an  infinite  integral.  A.  M.  Ostrowski.  Submitted  to 
the  Duke  Mathematical  Journal. 

(25)  Note  on  Vincent's  theorem.  A.  M.  Ostrowski.  Submitted  to  the 
Annals  of  Mathematics. 

(26)  On  two  problems  in  abstract  algebra  connected  with  Horner's 
rule.  A.  M.  Ostrowski.  Submitted  to  the  American  Mathematical 
Monthly . 

(27)  On  a discontinuous  analogue  of  Theodorsen's  and  Garrick's 
method.  A.  M.  Ostrowski.  To  be  included  in  the  "Proceedings 
of  a symposium  on  the  construction  and  applications  of  con- 
formal mapping,"  to  be  published  by  the  National  Bureau  of 
Standards  . 

(28)  On  the  convergence  of  Theodorsen's  and  Garrick's  method. 

A.  M.  Ostrowski.  To  be  included  in  the  "Proceedings  of  a 
symposium  on  the  construction  and  applications  of  conformal 
mapping,"  to  be  published  by  the  National  Bureau  of  Standards. 
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(29)  Formulas  for  complex  Cartesian  interpolation  of  higher  degree. 
H.  E.  Salzer.  Accepted  for  publication  in  the  Journal  of 
Mathematics  and  Physics. 

(30)  Polynomials  for  best  approximation  over  semi -inf inite  and 
infinite  intervals.,  H.  E.  Salzer.  Accepted  for  publication 
in  Mathematics  Magazine. 

(51)  Coefficients  for  polar  complex  interpolation.  H.  E.  Salzer. 

Accepted  for  publication  in  Journal  of  Mathematics  and  Physics 

(32)  Tables  of  zeros  and  weight  factors  of  the  first  15  Laguerre 
polynomials.  H.  E.  Salzer  and  Ruth  Zucker . Accepted  for 
publication  in  the  Bulletin  of  the  American  Mathematical  So- 
ciety . 

(33)  Formulas  for  numerical  differentiation  in  the  complex  plane. 

H.  E.  Salzer.  Accepted  for  publication  in  the  Philosophical 
Magazine . 

(34)  Formulas  for  numerical  integration  of  first  and  second  order 
differential  equations  in  the  complex  plane,  H.  E.  Salzer. 
Accepted  for  publication  in  the  Journal  of  Mathematics  and 
Physics . 

(35)  A bibliography  on  numerical  methods  in  conformal  mapping, 

¥.  Seidel.  To  be  included  in  the  "Proceedings  of  a symposium 
on  the  construction  and  applications  of  conformal  mapping," 
to  be  published  by  the  National  Bureau  of  Standards. 

(36)  Inequalities  concerning  ultraspherical  polynomials  and  Bessel 
functions.  0.  Szasz.  Submitted  to  the  Bulletin  of  the  Amer- 
ican Mathematical  Society. 

(37)  Summation  of  slowly  convergent  series  with  positive  terms. 

0.  Szasz.  Submitted  to  the  Journal  of  Mathematics  and  Physics, 

(38)  On  the  Gibbs  phenomenon  for  Euler  means.  0.  Szasz.  Submitted 
by  request  for  publication  in  the  12th  Anniversary  volume  of 
the  Acta  Scientiarum  Mathematicarum  (Hungary) . 

(39)  A remark  concerning  the  characteristic  roots  of  the  finite 
segments  of  the  Hilbert  matrix.  0.  Taussky-Todd . Accepted 
for  publication  in  the  Oxford  Quarterly  Journal. 

(40)  A recurring  theorem  on  determinants.  0.  Taussky-Todd.  Ac- 
cepted for  publication  in  the  American  Mathematical  Monthly. 

(41)  A problem  on  arc  tangent  relations.  J.  Todd.  Accepted  for 
publication  in  the  American  Mathematical  Monthly. 

(42)  The  condition  of  certain  matrices.  John  Todd.  Accepted  for 
publication  in  the  Quarterly  Journal  of  Mechanics  and  Applied 
Mathematics . 

(43)  The  condition  of  a certain  matrix.  John  Todd.  Accepted  for 
publication  in  the  Proceedings  of  the  Cambridge  Philosophical 
Society . • 
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(44)  On  conformal  mapping  of  variable  regions.  3.  E.  Warschawski . 

To  be  included  in  the  "Proceedings  of  a symposium  on  the  con- 
struction and  applications  of  conformal  mapping,"  to  be  pub- 
lished by  the  National  Bureau  of  Standards. 

(45)  Statistics  in  analytical  chemistry.  W.  J.  Youden.  Accepted 
for  publication  in  Transactions  of  the  New  York  Academy  of 
Sciences  . 

2.4  Reviews  and  Notes 

Symposia  on  conformal  mapping  and  Monte  Carlo  method. 

J.  H.  Curtiss.  To  appear  in  the  October  7th  issue  of  Science. 

At  the  request  of  the  following  journals,  staff  members  pre- 
pared abstract  reviews  of  technical  books  and  articles,  as 
follows : 

Mathematical  Reviews 

E.  F.  Beckenbach  - 2 

U.  S.  Quarterly  Book  List 
E.  ¥.  Cannon  - 2 

2.5  Miscellaneous  Publications 

(1)  Some  recent  trends  in  applied  mathematics.  J.  H.  Curtiss. 
Accepted  for  publication  in  the  American  Scientist. 

(2)  The  role  of  a statistical  consultant  in  a research  organiza- 
tion. C.  Eisenhart.  To  appear  in  the  Proceedings  of  the 
International  Statistical  Conferences. 

Operational  aspects  of  instrument  design.  Churchill  Eisenhart. 
To  appear  in  the  October  7th  issue  of  Science. 


(3) 


EXPLANATION  OF  PROJECT  DESCRIPTIONS 


The  project  descriptions  appearing  in  this  report  are  reproduced  from  the 
Project  Forms  used  in  the  project  control  system  of  the  National  Applied 
Mathematics  Laboratories.  With  a view  toward  making  this  report  more  useful, 
an  explanation  of  certain  terms  used  in  the  Project  Forms  is  given  here. 

Date  of  Authorization.  This  is  the  date  on  which  work  on  the  project  was 
authorized  by  the  Chief  of  the  National  Applied  Mathematics  Laboratories. 

Status.  Here  is  given  the  narrative  of  the  progress  to  date  on  the  pro- 
ject. Certain  descriptive  terms  are  used  to  indicate  at  a glance  the  nature 
of  the  activity  on  the  project  during  the  period  to  which  the  entry  applies. 
These  terms,  with  their  explanations,  are  as  follows: 

"NEW"  means  that  the  Laboratories  made  a commitment  within  the  quarter 
covered  by  the  report  to  work  on  the  project. 

"CONTINUED"  means  that  the  work  was  initiated  prior  to  the  quarter 
covered  by  the  report,  and  was  in  progress  during  the  quarter. 

"INACTIVE"  means  that  the  Laboratories  made  a commitment  prior  to  the 
quarter  covered  by  the  report  to  work  on  the  project,  but  no  work  of  any 
consequence  was  performed  on  the  project  during  the  quarter. 

"COMPLETED"  means  that  all  the  technical  work,  including  the  preparation 
of  manuscripts  of  the  final  reports  (if  any)  has  been  completed.  In  the  case 
of  tables  for  which  the  galley  proof  or  page  proof  is  to  undergo  extensive 
mathematical  checks,  the  designation  "Completed"  is  employed  only  after  these 
checks  have  been  performed. 

"TERMINATED"  means  that,  although  all  aspects  of  the  objective  were  not 
achieved,  it  was  necessary  to  terminate  the  project  due  to  circumstances 
beyond  the  control  of  the  Laboratories. 

Publication.  This  entry,  when  it  appears,  gives  information  as  to  the 
availability,  or  expected  availability,  of  the  results  of  the  project.  "In 
Manuscript"  means  that  the  results  have  been  written  up  and  are  available  for 
reference  at  the  Laboratories,  and  furthermore  are  in  a form  suitable  for 
photo-offset  or  other  means  of  reproduction.  In  the  case  of  "Completed" 
projects  for  which  manuscripts  of  reports  are  in  the  process  of  publication, 
further  periodic  entries  are  not  made  under  Status  or  Publication  to  record 
the  successive  steps  of  the  publication  procedure,  such  as  the  reading  of 
galley  proofs,  etc. 


